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- CALCULATION OF LOAD DISTRIBUTION IN 
STIFFENED .CYLINDRICAL SHELLS* 
By H. Ehner and H. Zoller 

Thin-wall e& shells with strong longitudinal and trans- 
verse stiffening (for example, stressed-skin fuselages and 
wings) may, under certain simplifying assumptions, he 
treated as static systems with finite redundancies. In 
this report the underlying basis for this method of treat-* 
mont of tho problem is presented and a computation proco- 
dure for stiffened cylindrical shells with curved sheet 
panels indicated. A detailed discussion of tho force dis- 
tribution duo to applied concontrated forces is givon , .and 
tho discussion illustratod by numerical examples which . re- 
fer to an experimentally detorminod circular cylindrical 
shell* Tho tost results reported by E. Schapits and C 
Krumling in a simultaneously appearing paper (reference l), 
confirm the reliability of tho method proposed. 



I. INTRODUCTION 



Tho determination of tho stress distribution in 
stro88ed-skln (or monocoque) wings and fuselages, accord- 
ing to tho elementary Navier theory for the bending of 
beams, or the Bredt theory for the torsion of thin-walled 
hollow bodies, is based on the preliminary assumption that 
the applied forces correspond to the elementary stress 
distribution. Large disturbances, which as characteristic 
additional stresses are superimposed on the elementary 
stresses, arise on the application of concentrated longin> 
tudinal forces - such, for example, as occur at a closed 
shell whore the continuity of the structure is interrupt- 
ed by a cutaway portion. (See fig. la.) 

Further deviations from the simple state, of stress in- 
dicated by the abovo theories arise on torsion and bonding 



"Zur 3ereehnung dee Kraft yorlaufes in verstelften Zylin- 
derschalen." Luf t fahrtf or schung , vol. 14, no. IS, 
Locembor 20, 1937, pp, 607-626. 
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by transverse forces if the loads are applied at the inter- 
mediate transverse walls which are restrained from defor- 
mation, or if the dimensions are not uniform along the ■ 
axis, or if the shell is rigidly attached at the end (fig. 
lb). Pinally, a small transverse stiffness of the bulk- 
heads (as, for example, in monoco-que fuselages) may, on 
tho application of arbitrary transverse forces and moments, 
affect the elementary stress distribution.. 

The computation of the additional stresses due' to ret* 
straint against deformation of a transverse section in 
torsion, has already been treated by several authors (ref- 
erence 2). A method of determining the additional stress- 
es arising. on the application of longitudinal . forces in 
long, thin-walled, cylindrical shells has been given in a 
papor by H. Wagner, and H. Simon (reference 3). For shells 
with strong longitudinal and transverse stiffening, a sim- 
plified shell model with redundant axial groups of forces 
has been proposed in a previous work, in which the force 
distribution is determined by a statically indeterminate 
computation.* This method has beon applied to the treat- 
ment of the special case of tho 8-stringor cylindrical 
sholl with double cross-soctional symmetry undor axial 
load (roforonco 5), In tho present paper the underlying 
principlos aro prosontod in a comprohensive manner and tho 
computation is extended to cylindrical shells of arbitrary 
cross section under bending, twisting, and axial loading. 
Tor the latter case there is investigated, with the aid of 
numerical examples, the effect of various stiff eners on 
the force distribution.. The method is further applied to 
a cylindrical shell whose stress distribution has beon ox- 

erlmontally determined for the caso of axial loading. 

Seo roforonco 1.) A corresponding computation for the 
caso of stiffened flat disks will appear in a later work 
in connoctlon with similar test results. 



.II. UNDERLYING PEINCIPLES 0! THE METHOD 
1. Simplified Modol of the Sholl Structure 



Tho cylindrical shells consist essentially of the 
shoot or skin, the longitudinal stiffeners (stringers), 



Tho same simplifications for the computation of stiffened 
shells are given by 0. S. Heck in his .paper (reference 4). 
His procedure, illustrated by examples, differs essentially 
from the one here presontod in the choico of the static 
redundancies* 
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and tho transverBO etiff oners (bulkheads ) . The longitudl** 
nal attf f eners are .assumed. -t,Q ho Inside -and attaohed to 
tho -skin, while the transverse stiffdners, assumed as par- 
allel, may olther ho attaohed to the shoot (for example, 
in the neighborhood of the positions of load application) 
or attaohed only to the inner aides of the stringers. Tho 
hxtVkheads may he in tho form of frameworks or flexurally 
rigid solid rings. The shell oross section is assumed to 
he simply connected, although the method may he extended 
to. multiply connected systems. lor the sako of simplic- 
ity, we shall restrict ourselves to cylindrical shells 
since tho disturbances of tho elementary stress condition 
in monocoguo fuselages and wings almost always extend over 
individual. small portions, which may approximately he con- 
sidered as of cylindrical shape. 

The transverse stiffening walls aro denoted in the 

positive x direction by 0, 1 k, .... n (fig. 2). 

They divide the system into bays, each of which is denoted 
by the stiff ener number on its. right. The stringers are 
numbered going around the positive direction from some ini- 
tial position 1, 2, .... J, m, the peripheral por- 
tions lying between tho otiff oners (shont or bulkhead por- 
tions) boing donotod by the higher stringer number as sub- 
script. All gpomotrical and structural magnitudes in tho 
sholl rocolvo a double subscript J,k corresponding to 
thoir position.- If tho transvorso section possesses axes 
of symmetry, then symmetrically lying magnitudes are given 
tho same subscripts. 

tfor the computation of the shells under consideration, 
tho following simplified model is used as a basis. The 
bulkhoads tako up forces that lie in the bulkhead planes 
only. The stringers aro assumed to be hinge-conne oted to 
tho bulkheads. In the case of stringers that go through 
and possess great floxural stiffness of their own, the ad- 
ditional stroBses may bo taken into account, as in tho 
case of frameworks, by the introduction of nodal -point" mo- 
monts as additional statio redundancies. In the skin it- 
solf a pure membrane- type of stress distribution is assumed. 

The most • important simplification is in conneotion- 
with the stress' distribution in the^skln and stiffenars. . . 
The transverse stiff eners are, for the moment, assumed to 
be attached to the skin; a generalisation will be consid- 
ered under section III, 4. Normal stresses will then be 
transmitted by the ' stiff eners only, so that the sheet pan- 
els act as pure tension fields. The normal stiffness of. 
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the skin is taken Into account by the addition of an ef- 
fective width to the stringer sections. Under these as- 
sumptions, in each of the sheet panels there is a con- 
stant "shear flow" 2 (t = s t) to which correspond linear- 
ly distributed normal forces in the longitudinal stiff en- 
ers within the panels, while the transverse stiff onor 
walls are stressed "by tangential peripheral loads, namely, 
differencos of tho shear flows in the longitudinal direc- 
tion. This state of stross agrees with the actual state - 
the more closoly the stronger the stiffeners as compared 
with tho skin. In the neighborhood of points of attach** 
ment and cut-out portions, these are the relations which 
generally exist. Tho reliability of the assumptions made, 
which had previously "been applied to the computation of 
statically determinate systems (bending beams) is con- 
firmed by tests. The assumption of constancy of the shear 
in each of the panels in the longitudinal direction may, 
in the case of curved shells, also be justified by the 
fact that generally the flexural stiffness of the skin, 
and stringers is small compared to that of the transverse 
stiffness of the bulkheads, so that within no panel can 
important peripheral stresses arido .in tho skin* 

This simplified shell model corresponds to a lattice- 
work, with redundant transverse frames in which the diago- 
nals are replaced by shoot panels under shear. For (n - 1 
intermediate transverse frames, each of which 1b q-fold 
statically indeterminate, the entire system with m string 
ors is [ (n - l) (n - 3) + (n + l)q_]-fold statically inde- 
terminate, and there arc further additional redundancies 
at points of attachmont. The statically indetorminato 
computation of the sumplifiod systom glvos the force dis- 
tribution "in the largo" and tho stross distribution in 
the actual shell design may then bo estimatod, taking into 
account the constructional dotails. (Soo tho numorical 
oxamplo in section 71, 2.) 



2. Principal Panel Systom 
Longitudinal Forces as Static Redundancies 

Tho computation of the statically indoterminate uodol 
shell procoods along linos similar to those developed by 
the first author in previous papers on space frameworks 



2 If 

The torm "Sohubf lusse , " translated as "shoar flow," is 
usod consistently in Gorman papers to denote . shear stross 
timos skin thickness. (Translator) 
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(roforonoo 6) and "box beams (roferonce 7). A suitable 
-principal system is obtained if each, intermediate bulkhead, 
in tho case of end restraint also at the end bulkhead, 
(m — 3) longitudinal attachments of the m stringers are 
loosened "by the insertion, for . example, of hinges. On ac- 
count of tho static redundancies in the transverse frames, 
this main system is still statically indeterminate although 
only stresses of the same transverse stiffenere may he su- 
perimposed. 

She redundant longitudinal forces in the resolved 
(m ~ 3) attachments of the stringers at a transverse 
stiffonor produce in the two neighboring hays (m - 3) 
stress distributions independent of each other. Such in- - 
dependent stress conditions may also be obtained in a 
"principal panel system 11 in the following manner.* By in- 
serting similar transverse stiffenere the bays are com- 
pletely separated and between- them there is Introduced at 
the stringers (m - 3) systems of longitudinal forces 
under equilibrium - systems which we shall denote briefly 
by the term "characteristic force groups" (Eigenkraf t- 
gruppen) (fig, 3). These systems then produce in the 
neighboring bays (m - -3) independent stress distributions' 
if none of the (m - 3) force groups may be represented 
by a linear combination of the remaining groups, tte then 
say that such groups are linearly independent. At each 
transvoree wall there are (m - 3) linearly independent 
characteristic force groups - for example, those consist- 
ing of a concentrated force and its "reaction forces" at 
three definite stringers not lying in a plane. The linear 
independence of (m - 3) arbitrarily set-up characteris- 
tic force groups may be established by the following cri- 
terion: There 1b formed the determinant of the (m - 3 ; ) 
order in which each row contains the forces composing a 
group at (n - 3) stringers where the remaining three 
stringers do not lie in a plane. The columns then con- 
tain the individual forces of the various force groups at 
some definite stringer. The (m - 3) characteristic force 
groups will then be linearly independent if - and only if - 
this determinant ' is different from zero. 

Example.* Figure 3 shows' three characteristic force 
groups 'f ot a 6-strrnger sheila- The- determinant , for ex- 
ample, of ' the Individual forces at the three upper string- 
ers will' be - 



See detailed presentation cited under references 5 and 6 
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- 1 0 + 1 

+ 1 -2 +1 

+ -1° - 1 + -1° 

2 7i 2y x 



and on account of ia different from zero, so that 

yi 

the force groups are linearly independent. 

Such linearly independent force groups at the k 
transverse stiffener and denoted "by ^i 1 k'« ^2fk» * " •^* t k» 
• ••» 3 ^3-3;k < are introduced as static redundancies and 
chosen as uniform in the axial direction - i.e., the com- 
ponent forces p j • of tiie unit state of ^ ^ are the 

same for all k's. In this representation of a principal 
system, it is to he noted that the double transverse walls 
are identical, so that the stress states in the right hulk- 
head of the k tn hay and in the left bulkhead of the 
(k+l)*^ . -fogy nus t -b e superimposed. In the system so chosen 
the stress states of the redundancies at the k* n trans— 
verse stlffener are superimposed only by those at the 
(k-l) th and (k+l) th transverse stiff eners at the k tn 
and (k+l} tJl hays, respectively, and by those at the 

(k-2) th - and (k+2) th transverse stiffener at the (k-l) th 
/ \ th 

and (k+1; bays. The redundancies at all other bulk— 

J.T. 

heads do not directly affect the redundancies at the k 
bulkhead« 

From an arbitrary system of (m — ?) linearly inde- 
pendent forco groups . X^, it Is possiblo to pass by a 
transformation. 

X H = °JA,1 Z l + °n,2 Z 2 + + c H,m-3 Z m-3 

(n = 1, 2 m - 3) (1) 

to new linearly Independent force groups Z^, where the 
determinant formod from the coefficients c [x,"o must be 
dlfferont from zero. This arbitrariness in the choice of 
the static redundancies may be utilized to obtain consid- 
erable simplification In the computation. For arbitrary 
force groups the statically indeterminate computation 
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loads, in the case of (n - l) intermediate transverse 
walls, to a system of (n - l) (m — 3) and n (m - 3) 
elasticity equations for froe and restrained end rings, re- 
spectively - the solution of which equations becoming vory 
tedious whon n and m- are large. An attempt is there- 
fore nr.de to choose the (m - 3) force groups in' such a 
manner that only the uniform groups in the axial direction 
affect each other. She elasticity equations then "break up 
into (n - 3) independent partial systems, each with 
(n - 1) and n equations, respectively, which with the 
principal system chosen, are composed of five members each, 
and in tho case of transverse stiffeners rigid in their 
planes, are composed of three members. 

The introduction of such "orthogonal"^ characteristic 
force groups is, in the case of an arbitrary number of 
stringers, possible only for a special shell form, namely, 
whero there is cyclical symmetry with regard to the geome- 
try and the elastic properties of the shell. The component 
forces Pj^' of each of the groups are then to "be set 

proportional to the ordinates of the sine and cosine curves 
with 2, 3, ••• waves over the shell circumference. & For 
m = 12 these force groups are shown in figure 4. (She 
circular functions with higher wave numbers lead to the 
same .groups .) With the aid of the orthogonality relations 
of the circular functions, it may easily be shown that tho 
groups of forces thus formed constitute systems in equi- 
librium, and that in their mixed displacement coefficients 
(boo section III, 2) tho individual contributions from 
stringers, sheet and transverse stiffeners vanish. They 
are thorofore orthogonal, as nay "be seen more simply from 
a consideration of the propertios of symmetry, since two 
forco groups - of which one is symmetrical with respect to 
an axis of symmetry of the system and the other antisym- 
motrlcal - do not influenoe one another. 



^The mutual effect of two such force groups is measured by 
the so-called "mixed displacement coefficient 11 (see III, 2; 
which is essentially a sum or integration of products of 
the corresponding stress values. She vanishing of such 
products is spoken of in mathematics as expressing tho 
"orthogonality," for example, of the vectors or functions 
cgncornod. 

Shis procoduro corresponds to the doformation method 
proposod by Southwoll for the computation of cyclically 
symmetrical space frameworks (reference 8). 
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The redundant, force groups may be eo ordered that 
with increasing subscript H more individual groups which 
aro in equilibrium over a smaller and smaller region of 
the circumference may he formed from their component 
forces; for example , in figure 4, for, the sine and cosine 
groups for \x = 2, equilibrium is possihle over the en*- 
tire circumference only, whereas the groups for \i = 3 
may each ho split up into two suhgroups which may ho in 
equilibrium over only half the clrcumf erenco. According 
to thq principle of St. Vonant, tho regions of influence 
of tho groups must decrease with increasing "order" of 
subscript [x, so that after a certain order thoy may be 
neglected. 

In tho case of a shell whose structure deviates only 
slightly from that corresponding to cyclical symne try, 
force groups are formed corresponding to complete cyclical 
symmetry and the slight mutual effects are neglected. For 
shells without cyclical symmetry, with only a few string- 
ers and" with sufficient axial symmetry of cross section, 
it is still possible to set up orthogonal forco groups* 
An oxample is that shown in figure 3 for a 6-stringer 
shell with cross-sectional symmetry with respect to the y 
and z axoB. Also for 12-stringer, doubly symmetrical 
shells, characteristic forco groups that have a negligiblo 
effoct on oach other may, as shown in section V, be set up 
for tho most important loading conditions. 

It is shown in the Appendix that for an arbitrary 
sholl construction, complete orthogonality of tho force 
groups can be attained in general only if the effect of 
the deformations of one of the threo shell- structure ele- 
ments is neglected - if, for oxample, rigid sheet or rigid 
transvorse stiffoner is assumed. The forco groups may in 
this case be set up, under the assumption of oqual panels, 
with tho aid of "principal axes transformation. 11 Shis 
proceduro, however, possesses only theoretical value. For 
most practical oasos it is possible to uso the simplo ap- 
proximate system considered in section T« An alternative 
is to restrict the computation to a few panels in the 
neighborhood of points of disturbance and solve the elas- 
ticity equations for arbitrary force groups according to 
tho iteration method described in section III, 3. 
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3. Decomposition of External Loading 
Principal Substitute Systom 

. Wo consider the loading cases In which there is a 
considorahlo disturhanco from the elementary state of 
stress. This will he the case whore concentrated forces 
are applied and also tho oaso of twisting and transverso 
honding of restrained shells or shells with "stepped di- 
mensions," and also the case where load is applied at tho 
intermediate "bulkheads. 

She loading of the shell hy concentrated axial forces 
1b generally reduced to the application of charact eristic 
force groups. lor this purpose there is subtracted from 
the external loading the forces corresponding to the ele- 
mentary theory and there 1b thus ohtained at the trans- 
verse stiffoner at which the load is applied, a force 
group which gives rise in tho shell to a state of stresB 
expressing tho differonce hotweon the olemontary and ac- 
tual stress. Bach force group lo acting at a transverse 
stiffonor nay, however, he split up into groups of tho 
form of the introduced redundancies , ?2» •••» ^n-3 

Xq = c;l Ix + °2 z 2 + •••• + °m-3 ^-3 ^2) 

since for tho component forcos of the groups at (m - 3} 
stringers, this equation represents a system of (m - 3) 
llnoar equations with the (m - 3) values , o 2 , 

c m-3 ' as unknowns. This system has a unique solution 
since the determinant formed out of the (m - 3) component 
forces of tho linearly Independent redundancies, according 
to tho criterion of section II, 2, is different from zero. 
The manner of this decomposition is indicated in figure 5 
for tho caso of a heading force group composed of four 
concontrated forces at a 12-stringer shell* The charac- 
teristic foroo groups and 7g correspond, to the co- 
Bine groups for p. = 3 and n = B in figure 4, and are 
considered more in detail in section V. 

-Tor tho external loading, or a portion of it, anothor 
principal system than that used for the static redundan- 
cies may he used in a statically indeterminate computation. 
This "suhstituto principal .system 11 is choson so that its 
state of stress corresponds as closoly as possihle to tho 
final one - tho statically indeterminate computation then 
constituting only an added computation; or a suhstituto' 
principal systom, for which tho computation of the load 
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coofficionts is as sinplo as possible, is taken. Both 
points of view, however, may not in general, he satisfied 
simultaneously. In applying a "bending force group consist- 
ing of four concentrated foroos (fig* 5 ) on a restrained 
shell, the four stringers at which the load is applied 
may, for oxample, serve as the suhstituto principal sys- 
ton and the load coefficients are then determined only hy 
integration ovor the nornal stress distributions in thoso 
stringors. The static rodundancios , however, then assuno 
largo values over the entire shell length. If, howevor, 
the bonding force group is reduced, as described, to tho 
application of characteristic force groups, the static re- 
dundancies die down to small values in the longitudinal 
direction. Then the shell is loaded by twisting forces, 
tho 'SJredt tube," and by transverse bending forces, the 
"ideal beam," respectively, are chosen as the substitute 
principal systems.. Tho statically indeterminate computa- 
tion then gives the deviations duo to the restraint against 
deformation, from the elementary stress condition. 

Shore occur, besides, under these loading conditions, 
disturbances due to the application at the transverse 
stiffeners, of forces not corresponding to the elementary 
theories. By subtracting the equivalent elementary stress 
distribution in the transverse stiff ener from that actually 
appliod, there is obtained an equilibrium system of extor- 
nal forces which deform the transverse stiffeners and thus 
affect the axial forco distribution. As a rule, however, 
the statically indoterninate conputation can be dispensod 
with if tho bulkhoads at which the loads are appliod are 
designed sufficiently stiff. 



III. PROOE73UEJB FOB THE STATICALLY IHDETEHMINATE OOMPUTATIOH 

1. Stress Distribution 



The following notation is usod for the dimensions of 
the cylindrical shell under consideration (fig. 2): 

ajj., length of k bay (bulkhead spacing) 

bj, width of sheet panel between (j - l)* 11 and J* 11 
*' strlngei. 
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m 

s J=l ^ J ' circumference of shell. 



m 



£ = E Sj t cross-sectional area enclosed "by shell. 

Tj fc, cross section of the J"* 11 stringer in the k th 
* hay (stringer section plus effective skin section), 

s« jg, skin thioknesB in k tn panel "between (j - l) th 
f and j ta stringer. 

For floxurally stiff ring "bulkheads with constant cross 
section along the circumference: 

■^S.ki cross Bection of k tn hulkhead. 

J S,kt moment of inertia of k** 1 hulkhead ring (with ro- 
spoct to radial deflection). 

The componont forces at the J tn stringer for the 

unit state of the redundancies *n,k are si'U-al for all 
values of k. The remaining symhols are sufficiently 
cloar from the toxt and illustrations. 

Ab a result of ^a t }c = It tne kth ( k + l )* h 

panels in the princi pal 'system are strossed. According to 
the simplified model sholl v normal forces arise in tho 
stringers that decrease linoarly to zero (fig. 6): 

Jik a k r J • ^J.k+l " ■ ak+1 r i K y 

(In each panel x is measured from the left hulkhead) and 
from 

sir** - <Vi,* 

there is ohtained for the constant shear flow in tho sheet 
fields .. . 

■ 

Tho shoar times thickness (shear flow) is chosen positivo 

in the direction of positive shoar stroBs, hence opposite 

that in tho work cited in roferonoo 2. in ' tho figures the 

roactions of tho shear flows on tho stiffonorB aro shown.- 
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(4) 



In particular cases 



may be determined 



from considerations of symmotry of the system and of the 
force groups; in general, however, there is required the 
additional condition that the shear flows at the trans-* 
verso stiffeners form a system in equilibrium whose moment 
therefore vanishes: 



where |. is the sector of the cross-sectional area from 

a suitable origin (for example, center of gravity of shell 
cross section) to the circumferential portion bj. 

Having obtained the shear flows, the stress distribu- 
tion in the transverse stiffener walls may be determined 
by considering the reactions of the shear flows as tangen- 
tial peripheral loads, the stresses in the (k - l) th and 

(k + l) tn bulkheads being opposite to the stress of the 

k th bulkhead from the shears of the k tl1 and (k + l) th 
panels, respectively. (See fig. 6.) Tho stress distribu- 
tion depends on the design of the transverse stiffener and, 
with static indeterninancy , required an additional static- 
ally indeterminate computation. Framework and solid con- 
nections are generally stiff enough, so that the changes 
in shape in their planes are small compared to the defor- 
mations of the stringers and sheet. The axial force dis- 
tribution will then be only slightly affected by them, so 
that for the statically indeterminate computation they may 
be assumed as rigid. Tho flexurally rigid bulkhead rings 
in monocoque fuselages, however, arc subject to greater 
deformations which in any particular loading condition 
must be taken into account. In goneral, they are threefold 
statically indeterminate, although important simplifica- 
tions result if axes of symmetry. of the cross section are 
present. The bulkhead ring computation for several shell 
shapes is given in section 7. 
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At the transverse . stiff ener 0, the applied charac- 
teristic force groups to which, the. axial force loading had 
"been reduced (section II, 3) stress only the first panel 
of the panel system. The state of stress 1b corresponding- 
ly given "by (3), (4), (5), where it is to be noted that 
the transverse stiffener 0 is loaded -only "by the shears 
of one T>anel. 



In pure torsional loading by moments M, 



bulkhead, the 
ment : 



.th 



at the K th 



panel- is aoted upon by a twisting mo- 



le- 1 

K=0 n - 



L 'k a vS„ «K ( g ) 
resulting, according to the Bredt theory, in a shear flow 



t<°> = - ^ 
k 2P 



(7) 



in the skin, constant ovor tho circumference and within the 
panel. The stringers are free from stress, the bulkhead 
k is loaded by the difference between' the actual applied 
twisting moment and the uniformly distributed elementary 
moment • 



For transverse force bending the determination of the 
"substitute state of stress" in the stringers and Bkin 
proceeds according to the well-known formulas of the "beam 
theory under the assumption of the simplified model shell 
(reference 2). In bending about the a axis, which 1b 
paBsed through the center of gravity of the shell cross 
section, the axial forces in the k* n ' panel are: 



and tho shear flows: 



s »i.k yi a 



. 3-1 

Ao) _ t (o) i.i *»* yi -s ,* 

'J.k 3 'l.k m 17 
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In the above equations B z is the tending moment, 
the transvorse force at the corresponding position, and 
(0) 

*1 k' when not obtained from the symmetry properties, is 
obtained in general from the condition: 

t (0) "h 

1=1 "i.k G 

for bending without torsion. In the transvorse stiffen- 
er nails there arises, as in the case of torsion, a stress 
distribution due to the manner of load application not 
corresponding to tho elementary theory. 



2. Strain Distribution 

Prom the stress distributions due to Z^ = 1, Z v = 

1, there are obtained in the usual way the strain o.r dis- 
placement coefficients 6„ v , that is, the virtual work 

which is done by the force group Z^ = 1 against the dis- 
placements due to Z v = 1 : 

6 1L „ = E / dx + Z b // - — 5-— du dx 

stringers^' Z 7 skin JJ s 8 G 



+ fi / 5 b — d +f i 1 — d (10) 

bulkheads L I J s B I Jg 5 J 

(du = circumferential element of skin, 
dug = circumferential element of neutral ring axis.) 

The bulkheads are assumed to be flexurally stiff rings 
in which, as in the case of the beam, the portion of the 
work done by the transverse forces, is neglected. For 
framework and solid connections, corresponding values are 
to be formed, the derivation of which will not be gone into 
since the changes in shape of such transverse walls are neg- 
ligible for the static computation. 



IT.A«C.A. Technical Memorandum No. 866 



16 



Te consider first the individual contributions .to the 
^jrir.tual work of the „ stringers t sheet,, and bulkheads sepa- 
rately. Por a stringer with cross ' seotion .P* and length 
a, the portion due to the nornal forces .1*^ n* 114 / £ 
which vary linearly "between the end values Pq , Pq ' 
and P<»*\ P< W > is: ' 

-jfe [pCi) (2P o ^) + P a -)) + p^) ( 2 P a -) + p(»)>] (ii) 

The displacement coefficient wf**' V ^ of the stringers of 
the k th hay due to X^ fc = 1 and X^ fc = 1 (or X^ = 

. 1 and x u,k-l = 



L,k 3 J=l E P 



A shoct of dimensions a. "b and shear stiffness s G con- 

(u) 

tributes the portion due to the constant shear flows t 
and t^: 

ah. t (ix) t (i>) ( 13 ) 

8 Gr 

so that the displacement coefficient w_ * of tho shoot 
portion of tho k tlx hay due to X^j. = 1 and X Vjk = 1 

< or X H,k-l = 1 * nd Xv.k-l 53 1 ^ is 5 



„<H.l>) b a k £ t<^ t< v > (14) 

B.k 12 j=l B j lr J« k 



The displacement coefficient w^» u ^ of the k th trans- 

S , k 

verso stiffener of each redundancy X^ ^ = 1 and X^ ^ = 1 

loadod only by the shears of the k* 11 hay in the case of 
floxurally stiff rings with constant cross section along 
the ciroumf orence , amount to 



xrCn.v) = O *(H,v) ( 15 ) 

3 . k 4a k * E J 3|k S,k 

with 
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whore , B^ U ^ denote tho tending nomenta ' 

the normal forces in tho k tn "bulkhead ring produced "by 
tho shears in tho k* n "bay due to X ji t k = 1 ani *v t k = 1 ' 

Tho coofficiont Kg^^ depends "besides on the typo and 

magnitude of the redundancies, only on the geometric shape 
of tho k th "bulkhead, the position of its neutral axis, 
and the stringer center of gravity. The coefficient var- 
ies little with the cross-sectional shapo and is given in 
soction 7 for several simple shell shapes. If the k tn 

+ "h 

"bulkhead is loaded only "by the shears of the k 611 and the 
(k + l) th "bays or only "by the shears of tho (k + l) th 
"bay, the coefficients aro, respectively l 

J*. AM) ^4 ffL K ( V ) 

a k+i 3 » k 4+i 8,k 

since the shears in the two "bays are proportional to ~- 

1- * 
and . 

a k+l 

From these soparate portions tho total displaconont 
coefficients 6^^ U ^ of the redundancies at the k tJl 
transverse stiffener may "be expressed in the form: 



6 k,k = [ W L,k + W L.fcflJ + L W B,k + W B,k+lJ 

In the longitudinal direction the redundancies k 
th 

at the k "bulkhead are combined with those of the 
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(k l)** 1 and (k + l)** 1 bulkheads, respectively , 
""""(^■Pjk— 1 an ^ k+1^ through superposition of the stress 

distributions in the k th and (k + l) tb "bays to form tho 

(u.w) (u.v) 
coefficients 8 fc k-1 and 8 fc ^ 1 , .respectively. Since 

the force groups at the different transverse stiffeners 
are. chosen uniform there hold the pynnotrical relations: 

- Cm- . v ) (v f p.) 
°k,k-l " °k,k-l " °k-l,k 3 °k-l,k 

and thoso "nixed displacement coefficients" nay ho ex- 
pressed in .torns of tho w values I 

B k,k-1 = 2 W L,k " w B,k 

- [if* (x + ffcl) + ( X + _*JL.) W (H.V)1 (l71} ) 

L a k > a k ' S,k-1 ^ a k+l' s » k J 

Purthor, hy superposition of the atross distrihutions in 
the (k - l)^ n and (k + l) transvorso stiffonor thoro 

arise tho nixod displacement coefficients 6/*\* „ and 

k f k-2 

g(l*i'0) t respectively, for which corresponding symmetrical 
k p k+2 

relations hold and which nay he expressed in terns of tho 
wg values:- 

a (|*.»> m fk-l (|^») (l7c) 
- k,k-2 a k S,k~l 

In a sinilar nanner are to he determined the coeffi- 
cients of the redundancies "by superinposing their .unit 
states upon the stresses due to the external loading" in 
the hay or substitute principal system according to (10). 
By the loading of the shell at the end stiff ener 0 with 
the characteristic force groups I t 0 , Z a 0 , o 

of the forn of static redundancies, only the redundancies 
at tho first and second transvorse stiffeners are affected. 

The load coefficients 8, „ . n ^ of the redundancies 

i.o * a . o 



» > 



i* ^p, a aue to a characteristic force group Z v Q = 1 
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are therefore: 



8 (n.u) = I w (n.») _ w (n.») r w Ci*.u.) + (i + wC^.t)! 

1,0 2 L,l B,l L S.O X aa/ S,l J 



ai (n.v) 



(18) 



2,0 a a S,l 

whore is to "be formed from the stresses of the 

s , o 

end bullchead 0 by the ahoara of the flrat hay due to 

1^ j = 1 and Z v L = 1. an arbitrary force group X 0 

at the transverse wall 0, and which according to (2) is 
split up in the form: 

n-3 

X o = u ^ 1 °u z u,o 

the load coofficionts of the redundancies Z (i) t and X^ j 
amount to 

>> = D i 3 =„. a'"'"', ^ c v , 6 ( ^> (19) 

1,0 v=1 V i,o a ( o v=i v 8,o 



In pure torsion by the moments at the bulkheads 

and according to the elemontary theory applied diatribu- 
tlvoly, there follows according to (7; and. (13) for the 

(n) 

load coefficients 6. of the redundancies X , : 

k,o li,k 

6 (n)_..^ g *k *J QQ _ *k+l » »k+l * J (|Q 
k.o 2 ^ J=la j f k a J.k 2 ^ *' =1 8 j,k+l & J.k+1 

(20) 

If the torsional momonts aro not applied distributive- 
ly, according to the elementary theory, there are still to 
be added the displacomont coefficients due to tho strossos 
in tho tronsvorse walla. 

In tho same manner in tho case of bending by trans- 
verso forcos. the load coefficients may, with the aid of 
formulas (11), (13), and (16), bo built up from tho indi- 
vidual contributions of stringers , shoot, and bulkheads. 
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3. Elasticity Equations and Solution In the General Case 

After the conputation of the virtual work coefficients 
the magnitudes of the redundancies are determined from the 
elasticity equations which express the conditions that the 
deformations of neighboring hays at the common "bulkhead 
should agree. lor arbitrary characteristic force groups 
as redundancies these equations have the form: 

. u=i Vk.k-2 T»fk-2 * ° kf k-i ^v.k-l * ° k§k A v,k * 

+ 8 k,k+l X ».k*l + 5 k,k+2 Z V f k+2; + 8 fcf0 = 0 (21) 

k 53 If 2, n - 1 and n; u = 1,2, «... m - 3 

with the "boundary values: 

2 H,-1 = ^n.o = 0 and = X^^ = 0; 

aEd X ji,n+1 s Z H,n+2 B 0 

for a freely deformahle end stiffener and end stiffener re- 
strained against deformation, respectively. 

The system of equations therefore breaks up, as shown 
in figure 7, for n = 6 and m = 5, into partial systems 
of five merahers each and which are formed in the direction 
of the principal diagonal from the displacement coeffi- 
cients of uniform longitudinal force groups, while the re- 
maining partial systems are composed of the mixed displace- 
ment coefficients of nonuniform force groups. • If the force 
groups are orthogonal to one another, these "mixed" partial 
systems, vanish. The values of the uniform redundancies in 
the longitudinal direction are then ohtained from Independ- 
ent 5— mem'oer partial systems which only contain uniform 
redundancies (principal equations). Tor the solution of 
these 5-nemoer elasticity equations - which In the case of 
rigid transverse stiffener walls are 3-memher equations - 
there have "been developed In static -structure computations 
a numher of suitable methods, so that it is not necessary 
to go int-o the matter any further. (See among others, the 
works cited- in reference 9.) For regular systems the so- 
lutions may "bo given In finite form as has "been dono in 
section IT, for equal panels. 



A 
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If the redundant force groups are not orthogonal to 
one another, it is generally possible to obtain the con- 
dition that the mixed coefficients of nonuniform groups 
he small compared to the coefficients of the principal 
equations. lor more accurate computations these approxi- 
mate solutions are improved step by step by substituting 
the values in the mixed partial systems and thus obtaining 
modified load members, for which the principal equations 
are again solved. This iteration process is most conven- 
iently carried through in "single steps" - i.e., the now 
approximate values of a group aro used directly to inprovo 
tho noxt group. Tho process thon always convorgos for 
elasticity equations (roferencc 10). 

After tho determination of the redundancies, the fi- 
nal force distribution may immediately be obtained "by 
superimposing upon the load condition in the principal 
system the corresponding multiples of tho unit states of 
the redundancies. 

If Pj°l is the force due to tho external loading 
in the principal system, the force Pj jg in tho J* n 

stringer at the k* n hulkhoad • is : 

n-S 

P, . = p(°) + Z p(M-) x . (22) 
Jik J.k ul=i J Htk 



Por tho shear flow tj k in tho j tn panel of the 
bay, thore is obtained: i 



-th 



IS— 8 



where t^°^ is the value in tho principal system. Corre— 

spondingly there is found the stress of tho k hulk- 
head from the redundancies 1^ , ^j lf ] E i and Z |x,k+1* ^ ne 

bending moment, for example, in the k* n bulkhoad ring is: 

- B<°> + V 3 M L x + (\ + -^A I , 

- B k + ul^i B k [ X H,k-l + V, 1 + *n,k 



B k = B. 



a k 
a k+l 



*l*,kfl] (24) 
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where 2v 0 ' is the value in the principal system, 

"IT the value due to- the -shears ■ of the k "bay as 
of X^ |k = 1. 

4. Extension of the Shell Model 

If the transverse stlffeners are not attached to the 
skin hut only fastened to the inner edges of the' stringers* 
as is often the case with monocoque fuselages,, the shell 
model thus far considered cannot he applied without modifi- 
cation. A transfer of the shears in the form of tangen- 
tial loads to the bulkheads is then posslhle to a limited 
extent only, on account of the weak attachment. In order 
to take this effect into account, the limiting condition 
is considered where no tangential forces at all can he 
transmitted from the sheet to the transverse stlffeners, so 
that there is a steady transition of the shear at the 
stiff oners. A state of equilihrium in the case of variahlo 
longitudinal, stresses is posslhle only through the setting 
up in the skin of peripheral stresses <r U f which ' corre- 
spond to the variable longitudinal shear stresses. In the 
case of curvod sheet with snail handing stiffness of its 
own* the radial components of those peripheral stresses 
must he transmitted to tho stringors which, as a result of 
their flexural stiffness, retransmit these forces to the 
transverse stlffeners. For the sake of simplicity, lot 
this radial load of the transverse stlffeners he assumed 
as continuously distrihuted over the poriphery as is ap- 
proximately the case for closely spaced stringers. ■ 

tfo consider first tho loading of the k* 1 * bulkhead 
for arhitrary axial shear distributions (fig. 8a). For the 
peripheral force P u = b 0"^, there is ohtained from the 

Bp,, 3t 

equilihrium equation — 31 + — = 0: 

3u 3x 

u 

where u 0 is a position at which the peripheral stress 
vanishes. The radial component q per unit of area at a 
position with radius of curvature r is then (fig. 8h): 



B 00 

k ■ 

a result 
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-1 p = I A I* dttf 



*o 

At the transverse stiffonors themselves hinge connection 
is aasunod at the stringers. The radial force E du ovor 
an olonont . du which acts on the transvorso stlffencr as 
a "reaction force" may then ho dotoriainod from the moment 
conditions. For tho k tJl hulkhead, as a result of the 

shear in the k* h hay (fig. 8c, x is computed from the 
left hulkhead in each hoy), thore is ohtainod: 

a k a k 11 

O 0 u Q 

Interchanging tho order of integration and integrating hy 
parts with rospect to x, there is ohtained: 



*k,k - i^tJ [*k*t*rj tdx ] du ■ i f C*(k) - V3 4u 

*o 0 u o (25) 

whore ls tho fhear flow at tho k tlx hulkhead and 



a k J 



'k 



V = at 7 * 41 
o 

is tho moan value of the Bhoar flow in tho k** 1 hay. The 
k tlL hulkhoad is further loaded hy the shears of the 
(k + l) til hay. For tho radial force B^k+i at position 
u t thoro is ohtained according to a computation similar 
to tho ahovo : 



*o 0 



k,k+l " o k+1 r 

if C- *ffc) + *k + l- d * (26 J 
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The total radial load of the k th bulkhead per unit cir- 
- — cumf orential distance ' is- there-f ore : 

*k - *k f k + ^k.k+i = 5 Ctk+i - V3 to . < 27 > 

If no tangential loads can "be transmitted, tho radial, load 
therefore dopends only on tho moan integrated values of the 
shears in the neighboring hays. 

To compute the system, tho "simple shear field scheme u 
according to section II, 1, is used as a starting haoiB. 
The longitudinal stiffnesses of skin and stringers are 
combined and in place of the variable shear in the longi- 
tudinal direction, there is taken the mean integral value 
within each bay and the corresponding linear force distri- 
bution in the stringers. In the statically indeterminate 
computation for the redundant axial force groups, there 
are neglected at each transverse stiffener wall the dis- 
placement contributions of the skin due to the peripheral 
stresses, and their effect is taken Into account only 
through the varied loading of the transverse stiffenors 
with radial forces according to (27). There are then ob- 
tained othor displacement coefficients ^g^^ than for 

tho case of tangential loading. The difforenco depends, 
howovor, as will be shown by an ozample in section 7, 1, 
only on tho ex centric position of skin bulkheads. Sinco 
no abrupt discontinuity of the shear oan occur at the 
transverso walls, tho computed "stepped" curvo of tho moan 
shear flows in tho longitudinal direction must be . smoothed 
out by a ao ntinuous curvo in such a manner that the mean 
value within the bays remains tho same. (See fig. 8a.) 
Similarly the Btringor forces which are obtained by the 
statically Indeterminate computation: at the transverse 
stiffonor. walls must be joined by a corresponding continu- 
ous curvo. It is possible,- howovor, to consider tho stiff** 
noes of tho bulkheads as distributed over the shell longth 
through the flozural stiffness of the longitudinal stiff- 
enors and in tho solutions for tho redundancies pass to 
tho limit of bulkheads. spaced Infinitely close. Thoro is 
thon obtained directly for tho stringer forces a continu- 
ous axial distribution.. Jor tho case of loading of a vory 
long shell composed of equal panels by means of characteris- 
tic force groups at an end bulkhead, these solutions are 
given in IT, 2, in the table of formulas 3. From the shear 
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distribution the peripheral forces in the skin nay he de- 
termined as has he en done in the work of E. Schapltz and 
G-. Znmling, cited in reference 1, starting from the experi- 
mentally determined longitudinal stresses. 



17. "SIMPLIFICATIONS FOE THE CASE OF EQUAL PANELS 
l a Displacement Coefficients and Elasticity Equations 



Of the simplifications of the statically indeterminate 
computation with uniform dimensions over the shell length, 
the case of geometrically and -structurally equal panels will" 
he considered. Systems of stepped dimensions may he ap- 
proximated to the above case if the panel dimensions in the 
neighborhood of disturbance positions (cut-outs, deforma- 
tion restraints) are usod as a basis for tho computation* 

In tho notation the subscript k is dropped; the equa- 
tions for the stress distributions havo the same forms - 
the w values according to (12), (14), (15), from whidi 
the displacomont coefficients are formed being tho samo 
for ?,11 k'o. It is convenient for the computation to 

- . (H,v) _ (n,l>) (l*,v) . 

brook up w L and w B as well as w g in 

(15) into a coofficiont k(^« v ) and a factor u) which 
depends only on tho system dimensions. In addition, wo 
writo I 





U)_; w 



(H.v) _ _1_ K (n-.v) 
B ~ aE B 



U). 



B 




U) 



S 



(28) 



with 



id 





(29) 



Thon, according to (12) and (14): 
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4 J=l Sj- J j 



(30) 



(**» W ) J -U* J ,1 

Kg is olJtainod 



Tor floxurally stiff "bulkhead rings 
from (16) as: 

«S (P " V) = B^> B< V > aus + lg « ^ v (lO *(") dng J (30 ,) 

-( t j£ ^ and t( v ^ are the shear flows in the k tlx panel 
duo to k = 1 and * v ^ a 1, respectively; b' , 

(n) (u) (v) 

IT and B , N are tho "bending moment and nornal 

forcos in the k** 1 "bulkhead duo to these shoars.) 

Tho redundancies at the k tn bulkhead are similarly 
combinod with thoso lying to the right and to the loft: 

(H,l>) (n,v) (n,u) (H,v) 
6 =6 : 8 = ft 

k,k+l k.k-l* k,k+2 k,k-2 



According to (l7a,b,c), tho dlsplacomont ooofficionts using 
thp abovo notation aro : 



3D 6 (11 ' U) = BK<V" V > U) + 2K<^ V > U) + | K^« w ) 
k.k L L B B '2 s 



a S 8 



a X 8 



k,k-l 

(H.v) 
k,k-2 



S 



„ (»*|V) (n,v) <i*.v) 

L LB B S S 



1 (jifW) 

4 *S "8 



(31) 



Similarly simple expressions are obtained for the load co- 
efficients. In applying the force groups z p ( o " 1 of 
the form of redundancies at the "bulkhead 0, they may "be 
expressed in terms of the displacement coefficients of the 
redundancies since 
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(v-,v) (H.V) (p.,v) (H.,1?) . 

8 =8 + ; ■ 6 =6 (32) 

1,0 k,k-l k.k-2 2,0 k,k-2 

If the end bulkheads have the same elasticity as the inter* 
mediate stiff eners. If the end bulkheads 0 and n, on the 
contrary, are very stiff compared to the intermediate hulk- 
heads (rigid in the Uniting case), then the load coeffi- 
cients are: 

6 1,0 ■ S k,k-1 + 28 k,k-2 S 6 2,0 = 8 k.k-2 (33) 

In this case the following displacement coefficients of the 
redundancies also change I 

(H.v) (n.v) (tx.u) '(n.u) 

1,1 n-l,n-l k f k k,k-2 

The load coefficients for the momont coefficients 

at the bulkheads are according to (20), since t. , ,_ = 

J v Jc+1 

- tCi*) = ^ tjiO. 

5 II*; - _JE E _J J (35) 

k,0 21 J=l sj 0 

With hays of equal dimensions therefore, the load coeffi- 
cients of the redundancies vanish at the unloaded trans- 
vorso stiffener walls. The samo holds true for the bond- 
ing undor transverse forces. 

The elasticity equations (21 ) may under these simpli- 
fications be considered as a simultaneous system of (m - 3) 
linear difference equations of the fourth order with con- 
stant coefficients whose solutions can be obtained without 
too much computation work in simple cases. (See reference 
11.) In the case of orthogonal characteristic force groups 
(or force groups which affect one another to a negligible 
extent), this system breaks up into (m - 3) independent 
difference equations of the fourth order of the form: 

8 k,k-2 X k-2 + 8 k.k-l Z k-1 + 8 k,k Z k + 

+ 6 k,k+l + 8 k,k+2 *k+2 = - 8 k,0 (36 > 
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where the shorter notation- X^, 6 k fc etc., has "been writ- 
'-ten for- 2^ j6 r, a^^J*) etc. She solutions of these differ- 
ence equations must, on account of the incompleteness of 
tho first and last elasticity equations (soo fig. 7) satis- 
fy certain "boundary conditions. 

Tho obtaining of solutions in finite form is poBsihlo 
for a froo "bay system - i.e., for a shell with nonrestrainod 
end "bulkheads. In tho case of an end "bulkhead n re- 
strained against doformation, the rodundanoios at tho 
fixed ond are dotorminod, using tho froo hay systom as a 
statically indotorminato principal system. If in tho froo 
X o ) 

system are the redundancies due to the external 

(n) 

load, X^ the redundancies due to X^ = 1, the final so- 
lutions are of the formt 

From the last elasticity equation, there then follows: 

T - 6 n.n-2 *n°2 + 6 n,n-l + 8 n,o . 

a n f n-2 *n-2 + 6 n t n-l "Si-l + 8 n,n 

The magnitude of the longitudinal forces at the fixed 
end restrained against deformation (fig. lh) must "be de- 
termined, for example, hy a statically indeterminate com- 
putation. The redundancies X£° ' are generally negligibly 

small in the case of uniform structure of the closed shell. 
Tho rodundanoios X^** aro oomputod according to tho moth- 

od givon for tho application of longitudinal forces. The 
rodundanoios X n at tho rostrainod end can thon "be dotor- 
minod according to (37). Prom a knowlodgo of tho force 
distribution on the application of concentrated forcos, 
tho main deviations from tho. olomentary stross conditions 
can thoroforo ho dotorminod, and for this reason this load- 
ing condition will ho considerod in detail. 
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2 # Solutions in Pinite ffarm of tho Difference Equations 

To solve the olasticity oquations (36) considered as 
symmetrical difference equations for uniform force groups: 

X k-2 + 27 X k-1 + 2 P 4 + 2Y X k+1 + Z k+2 - " ^ < 35 > 

(k = 1»2, , n - l) 



with 2Y =, 2^ = JkJE-! H k = 'fc.O (39) 

6 k t k-2 8 k,k~2 8 k,k-2 

we start with the general solution of the homogeneous equa- 
tion (Tljj. s 0). An exponential substitution leads to a 
characteristic equation of the fourth degree. Using hyper- 
bolic and circular functions, respectively, there is ob- 
tained the following result: 

She general solution of the homogeneous equation is 
composed of four independent particular solutions whose 

form depends on tho value D = 3(P g iO .. The case occur- 

ring most in practice, namely, Y < 0 'is' assumed; for 
Y > 0, the general solutions are to ho multiplied "by 

(-D\ 

Tar D > 1, the solution of tho homogeneous equation 

is_: 

= Oi cosh k * cob k X + Oa cosh k \|/ sin k X 

+ 0 3 sinh k * coo k X + G 4 sinh k \|f sin k X (40 ) 

1 

Tho arguments V and X satisfy the bou ndary/ condlti-ons: 

fCf r J . . .* \ . 

cosh * cos X = I sinh ♦ sin X = j3y| /5 - 1 (41 ) 



whonco follows: 



\j/ = J cosh- 1 (A + B) A = P = 3 " » .O^ 

with 

X - § (cos- 1 ) (A - B) B = y(fl-j-l)" - Y 5 " (42) 
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(Of the multiple-valued inverse trigonometric functions, 
only the posi.ti.ye principal values are. taken.) 

If D < 1, then in place of the oircular functions, 
there occur the corresponding hyperoolic functions: 

z k ~ °i oo sh k >|/ cosh k p + O a cosh k \|r sinh k p 

+ O a sinh k i|f cosh k p + 0« sinh k \|f sinh k p (40 «) 

with 



I Y 

cosh y\f cosh p = I — 



; sinh >|r sinh P 3 j|| D < 41 1 ) 



* = | cosh- 1 (A + B); P = | cosh" 1 (A - B) (42 ' ) 

(A and B as in (42)) 

lor Dal (double root of tho characteristic equation) 
X and p, respectively, are equal to zero, and 

Y 



i|f = | cosh" 1 (A + 3) = cosh"" 1 



2 



In addition to cosh k ^ -and sinh k 'Jf, k cosh k ^ and 

k sinh k A> are also solutions, so that the general solu- 
tion of the homogeneous equation is: 

I k = Oi coah k * + Os k cosh k i|f + 0 3 sinh k * + 

u + 0 4 k sinh k * (40 - ) 

Tho complete solution of the differonce equation 1b 
then made up of tho solution of the homogeneous eqiiation 
with four arbitrary constants and an arbitrary particular 
solution of the nonhomogoneous equations In the case' of a 
system witfh unrestrained end, the following boundary con- 
ditions must, on account -of the incompleteness of hoth 
the first and last elasticity equations, he satisfied: 

V = 0; 2^,1 = 0; X n = 0; Zq+i = 0* (43) 

. . 4 

There are thus ohtained four linear equations for the con- 
stant's Ci t C 8 , Cg , C« in the general solution, and after 

these have "been computed the solution of the elasticity 
equations is ohtained in finite form. 
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As en example, ire consider the case of constant load 
coefficients Tl^ = T| = constant, as is the case- in loading 
the cylindrical shell with constant moments at all trans-* 
verso 8tlffoners and where thore is no restraint on the 
deformation of the ond sections. Further, let Y be as- 
sumod as < 0 and D > 1. /A particular solution of tho 
nonhomogeneous difference oration is: 



* = ~ 



2 + 4Y + R3 



to which must he added the solution (40 ). After determin- 
ing the constants C from the "boundary conditions (43) 
(this formal computation is omitted), there is obtained 
the solution: 

= V " [{b 1 ^ (k + l) * sin (n - k + 1) X 

- sinh k \|/ sin (n - k) x) . 
+ -[sin (k+1) X sinh (n-k+l) \|/ - sin k X elnh (n-k)\|/}1 
where II* = sinh (n+l) \|f sin X + sin (n+1) X sinh ^ 

(Por 33 < 1 the hyperbolic functions are to be substituted 
everywhere for the circular functions with the argument p 
in place, of X, while for D = 1, the functional symbol 
sin and the argument X are to be dropped.) 

Tho sotting up of a solution of this type is possible 
for further simple loading conditions and it is also possi- 
ble to tako into account modified dimensions of tho end 
bulkho?,do by corresponding end conditions. Since tho gon- 
eral formulas aro not very .explicit, however, it is bottor 
in any individual cr.so to use a numerical nothod, or in 
the caso of a fow oquations, to solvo diroctly by tho woll— 
known elimination or itoration method. Por tho limiting 
case of rigid bulkhoads, solution formulas havo"boon sot 
up for a sorlos of load conditions in the work cited under 
reference 6. 



In applying independent force groups of the form of 
static redundancies at the bulkhead 0, the load coeffi- 
cients dlffcront from zoro of tho first and second equa- 
tions may bo detormlnod by the boundary conditions and 
simple solutions in finite forar may be obtained. Accord- 
ing to formulas (32), (33), (."4), wo must have 
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Z 0 = 1; = 1; 2^ = 0; X n+1 = 0 (44) 

for tho caso of elastic end bulkheads, and 

I 0 = lj + =» 2; = 0; + J^ +1 = 0 (45) 

for the case of rigid ond bulkheads. 

By neons of those conditions, tho four arbitrary con- 
stants 0 in the general solution of the homogeneous 
equation: aro determined. Tho final results aro preeontod - 
in tho table of formulas 1. She solutions for all cases 
may be expressed in a gonoral form through auxiliary func- 
tions g fc and tho conditions Z Q = 1; = 0 may bo in- 

modlatoly verified. 

In tho table of formulas 1, there are given in addi- 
tion the solutions for the limiting case of rigid bulk- 
head iralls. She differonoe equation then roads: 

+ 2a I k + X^ =0 (46) 

with 2a = — * and tho end values X 0 = 1, = 0 . 

8 k f k-l 

Tor 1, the equation has the general solution 



X k = (±l) k [Oi cosh k q> + 0 a sinh k <p] 
(upper sign if a< 0) with the argument 

q> a cosh -1 | cc | (47) 
For J a =1 the solution is: 

Ik = (±D k [Ox + % k] 

There are thus obtained for tho above boundary values the 
solutions : 

_ f . , *k s inh (n - k) cp 
= I* 1 * - 8lnn n cp 

and (±l) k (l - i k) 
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Form of solution: X k - (± if ** " ^ ; "" H " r ^ ^ ' * ° ° r * <0 

9, -9A lower s/gn ifyX) or <*. »o 


Case 


Inter mediate butkheac 

elastic 


Is e/astic, end bu/h heads : 
rigid 


D > 1 
(A-B<1) 


g„ - sm(k*1)xsinfi(n+1-k)y 
+ sinkx sinh(n -k)y 


g K ' sin(k+1)x sinfi(n+1-k)y 

7 2 sin kx •sinh(n-k)Y+sin(k-1)xsinh(n-f-k)y 


D- 1 
(A-B-1) 


g„ - (k+1) sinh(n* 1-k)y 
7 k s/'nh(n-k)f ' 


gk*. (k+1) •Smh{n*1-k)y 

7 2 k 3iah (n-k)y* (k-1)s,nh(n-1-k)y 


D<1 
(A-B>1) 


g k is/nb(k+f)<? s'mh(n+1-k)y 
7 s'tnh k<> 5ir>h (ll-KJy 


g k -5inh{k+i)Q sinh(n+1-k)y 
r2sinhkp s'mh(n -k)%> +s<nh{k-1)p''nh[n-i-k)y 


Limiting case of rigid bulkheads :g h 'sjnh(n-k)<p fOrjocj* /; g H = ri-k fOr jocj • 1 



Table I.- Solutions of elasticity equations for applied characteristic 
force group Xo=~i for shells of finite lengths (end bulkhead 
freely def ormable) . 

i 



Case 


Infermed/afe bulkheads, 
elastic 


elastic, end bulkheads: 
rigid 


A-B < 1 


X k --(±lfe^[coskx* c -^s)nkx\ 




A-B - 1 


X k -(*lfe^+[l + (!+e-*)k] 




A-B > 1 






Limiting case of rigid but M heads (a> s - 0): X k -(±1) k e' kf 


Upper sian if x s u) s >2x L (0 L -x B to ei Lower sign if x s (0 s < 2 ' Xfi; L -X B co 5 


Constants. y - -l-sirn rnstfA +R) 

. _ 8x L u L +2 Xe v B +x s u> s i 
A . x s u s X • J' arc cos (A-B) 

A ■ — , 4 

IXbCOb * 4x l C*>l 1 (p - 4- ore nashloLl 
/0£/ Ix b cj b -2x L cotl \ 



Table II- Solutions of elasticity equations for applied characteristic 
force group Xq^ 1 for infinitely long shells. 
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Case 


Intermediate bulkl 
elastic 


ie-ads elastic, £>f?d bulkheads; 
rigid 


A-B<0 


X(x) - e'^ [cos xx + fsin%x] 


X(x) - e VJt [cosxx * {(%-j)sin xx] 


A-B-0 


X(x)-e- f '[l±yx] 


X(x)' e' f *[l +t?*J ' : i 


A-E>0 


X(X) - e"^ x [cosi^X * ys'mhyx] 


Xcw- e**£osh?x fjCf+fJ^fx] 


Limiting case of rigid bulkheads : X(X) * G * 


Constants x g CO B 
x s o5 s 


v - 1/a+b ' 
x - YB-a' 
F - 1/a-B ' 



tfafelelXI- Distribution of force group 2(0) = 1 for infinitely long 
shell with bulkheads lying infinitely close. 



CO ~ 

Ja/ues 




X - 

Va/ues 


Arbitrary doubly symmet- 
rica/ cross section. 


Circular cross section with 
uniform stringer Spacing 


V- 
li 

0) 

Q 


a 3 - 


> *8* ~J~ 
*s~ Astir** + W8*i t ffi+ j£) 

-5,sh7oM s ] ro' z 


V. 

V 

a 




- 103,0099 \A S ]- -TO' ^ 



TfcblelV.- Displacement factors for six-stringer shells of equal hays. 
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The solution formulas simplify considerably for the 
.limiting case of inf inite., shell. Jength i.e. , n — > co for 
constant finite bulkhead spacing a. The results for this 
limiting condition are presented in the table of formulas 
2. These solutions may he applied to shells of finite 
lengths .greater than their perimeters, without appreciable 
error, as is shown with the aid of examples in section. 71, 
1 below. The solution becomes particularly simple for the 
limiting caso of rigid bulkheads. If in a, 0, Y, the 6 
values are expressed in terms of the K and cju values 
corresponding to (31), these coefficients need not be fur- 
ther specially determined - the constant values of the 
solution formulas being computed directly from the K and 
U) values, as indicated in the table of formulas 2. The 
arguments ty,X are taken from tables of functions 7 with 
the natural numbers as argument. The condition 

3) -~- 1 and A - B -3- 1 

is equivalent to 

12 (K L U) l ) (Kg U)g) (2K L U) L - K B U) B ) a 

and the condition Y 0, according to which the sign 

in the solution formulas is determined becomes, in terns 
of the K and u> values: 

For some simple shell shapes the K and u> values are 
given in table of formulas 4. 

A further limiting case of importance is that of in- 
finitely close bulkhead spacing for which the total trans- 
verse stiffness of the system is assumed to remain un- 
changed. Ef the bulkheads, for example, are rings with 
the bending stiffness EJ g , the bonding stiffness EI g = 
= EJg per unit length in tho axial direction is to remain 

constant. The difference equation then goes over into a 
corresponding differential equation, the independent force 
group X]g becoming a continuous function 3C( Z ) of the 
axial coordinate x measured from the loading side. Again 

— - ; 

For example, the tables of circular and hyperbolio func- 
tion of I. Hayashi, or "Hutte," vol. I. 
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we omit the computational details of passing to the limit 
and present the final solutions in the table of formulas 3 

The case .corresponding to V > o drops out, the cases 

▲ B -~- 0 correspond to D 1, the arguments are de 

termined by square-root expressions, and in place of the 
K and u> values there enter corresponding Tf and in 
values which no longer contain the bulkhead spacing a. 

If the shear deformations are neglected (k b u>£ = 0), 
then \p = X, x and the solutions are: 

— ipx — _ . 

X (x) e 6 * cos * x + Bln * x ' 

for elastic end bulkhead, and 

= e - ^ cob T x 

for rigid end bulkhead, with 

«/ «s ro s 

They have the same form as the function f e (z) for the 
simple disturbance loads in the work of Wagner and Simon, 
mentionod in reference 3. The argument values likewise 
agree for . corresponding stress distribution, as may be 
shown by the example of a box beam of sides b and c and of 
oonstant wall thickness s. From the work cited under 
.reference 7 (pp. 76-77), there is obtained for equal pan- 
els: 

k - ga a it ui ■ b B c a (b + o) 
" L ' b (b ■* c) 1 ** ^ - 24a~Jg 



so that 



■- 6 _ ,„ b a o a (b + c) 

K T = ; rt Kr da = — f — = 

L B ( b + c ) s s 24 I s 



Ye thon have: 



7 a 8 /b ±-±=2 

V /b^ J all* 



13' 856 
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in agreemont with the corresponding U) e value , according 
to tahle 2 of reference 3. 

7. CONSIDERATION OF SEVERAL IMPORTANT SYSTEMS 



For several simple shell shapes with few stringers 
and whoso cross sections are symmetrical about two perpen- 
dicular oxob, there will "be computed the stresses and dis- 
placement values. The transverse stiffeners are taken to 
"be floxurally stiff rings of constant section over the. en- 
tire circumf oronce and with radii of gyration small com- 
pared to the diameters. Shells whoso cross sections devi- 
ato llttlo from the condition of douhlo symmetry (nono- 
coque fuoolagos) may "be treated in approximately the same 
manner as theso simple systems. 



The simplest statically indeterminate shape of shell 
is the 4— stringer system with intermediate transverse 
stiffeners. The stringers are assumed not to lie on the 
axis of symmetry. As redundancies, there are set up at 
each intermediate transverse stiffener k and at the end 
restraint, a force group Zjj. symmetrical with respect to 

the center and denoted "briefly as a "convexing force grouj 
The unit state = 1 is shown in figure 9, and the re- 

actions of the constant shear flows ti and t a on 

stringers and rings, which for the present are assumed to 
he attached to the skin, are indicated for the k. hay.. 

In the stringers according to .(3) are set up axial 
forces which drop off linearly from 1 to 0, and for the 
shear flows there is ohtained according to (4) and (5 ) I 



1. Four-Stringer Shells 





*k li + Is' 2 ' k ~ * k Zi + la 



> 



(48) 




; t 



2.k+l - a . 



k+1 



1 




The stresses of the flexurally stiff, threefold stat- 
ically indeterminate Bulkhead rings (On the computation of 
bulkhead rinss, see among others, reference 12) may in the 
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case hero considered of double symmetry and loading symmet- 
rical, about .the.„center , "be obtained without statically inde- 
terminate computations, since at the four dianotrol points 
on tho axos tho "bending moments and the -normal forces van- 
ish and tho transverse forcos Qy and Q z at these points 
are obtained from tho equilibrium conditions. Lot ub. con- 
sider tho k** 1 ring undor tangential loading by tho shears 
of the k th bay (fig. 9b). Substituting from (48), we havo : 

k 7 V (49) 

With the aid of the above equations, the distribution 
of the "bending moments, and tho normal and tho transverse 
forcos ovor the entire ring may roadily be obtaino.l. Tho 
bonding uonont B^, for oxamplo, at a point 7s,k' a S,k 
of the neutral axis of the ring soction is: 

B t = (L 2 g v . ti v (2f x ) in poriphoral range b x ] 
k *y,k 5,k l.k V (so ) 

B k = ** ^z.k y S,k + ( 2f »* in P° rl P ]lC!ral range b a J 

(Bending momont positiyo if the extreme outer ring fibers 
are under pressure.) fj_ and f a are tho hatched areas 
indicated in figure 9b. The ring sootion may bo dif fo ront 
for oach bulkhead, and hence also the coordinates ys t ic» 
z g k of the neutral axis. The coordinates y x , z a of *ke 

stringer center of gravity are assumed constant over the 
entire shell .length. .The ■ k* n bulkhead is'in'-a similar 

th 

manner loaded by the' shears of the . (k + l).. panels it is 
only necessary to substitute a^+i a k" 

If the shell cross section is an ellipse: y = r y cos cp, 
z = r z sin cp; y x = X L r y coo cp x ; z t = X L r z sin cpj , the 

coordinates of the stringor conters of gravity and the shapo 
of the noutral axis similar to that of - the. circumference 
around tho -skins y g k = ^g^ r y cos cp; B s t ~k = 
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*S,k r z flln *P 8 (where Xj,, Xs t k are absolute numbers less 
than 1), there is obtained, after a short computation for 
the bending moments t 



=k = ~ [(*L *S,k c °* *i) sin qp - (l - ^) q>] 

in peripheral range b x 

- £(*L x S,k sln *Px ) cos cp - - <p)J 

in peripheral range b a 



> (51) 



r v r 
B k = -J- 



For a circular cylinder, cp and cp 1 are the angles sub- 
tended at the center (fig. 10), and 

r s k r L 

r y = r a = r; X s k = * L = y 

In this case the normal forces Ujj. may readily be 
obtained? 



¥ t " a k x 



co 



s q>x sin cp for 0 < cp < cp x 



(52) 



ffk = t- *L sin cp, cos cp for cp < <p < 2 



If the rings are not attached to the skin the same 
"mean" shear stress distribution of the extended shell 
model. is assumed according to (48), while the rings are 
loaded by distributed radial forces. According to (27), 
section III, 4, this loading B k of the k*h ring per 
unit distance around the circumference and resulting from 
the circumferential forces in the k tn and (k + l) th 
bays due to X k = 1 , is: 



E k - 



8 

ThiB assumption holds true approximately for constant 
ring section, since the -neutral axis is displaced inward 
for strong curvature. (See reference 13.) 
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In the first and third quadrants, the radial load is there- 
fore compressive while in the second and fourth quadrants, 
it Is tohe'ilo. '2^ is split up into two parts proportional 



to i 



and — = — . 
a k a k*l 



The tending moments and normal forces 



in the k th bulkhead due to the contribution from the k th 
hay are then: 



a k 



- £(*;& x S,k «b <Pi) Bin ep - X S(k (l - qpj 

for 0 £ qp £ <P X 

B k°^ [<*L *S,k 8la ^5 cos <p - X Sfk ^ - cp)] 



> 



(53) 



TT 



for qpj < cp < - 



> 



(54) 



Hfc = ^ [(X t cos (pj sin qp - (l - ~i) qpj 

for 0 < qp < cp x 

5k = ^ [ ( *& 8in ^ 008 * " (I " *)] 

for qp x < V < § 



The bonding r»omsnts differ from the values in tho caso of 
tangential .loading according to. (51) only through the "ex- 



oentricity 1 



occurring as a factor of tho 



"S.k r 

■ second meribor, and they aro therefore larger since 
Xg | l c < 1» The normal forces, however, are considerably 

smaller, a-j nay bo seen ^by comparison with (52). 

The displacement coefficients, corresponding to (17a. 
b,c), are made up of the w values, defined in (12), (14;, 
(15) I- 



»L,k - ? 



l.k 



(55) 



r B 
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The value of - Kg k according to (16) for the case of 

cross sections of general shape, may he determined only "by 
numorical integration; - for example, hy the Simpson rule.. 
Tor tho circular "bulkhead the integrals may he explicitly 
evaluated. In the case of tangential shear load with the 
stressos according to (51) and (52), thore is ohtained: 

g 

«B.k - V.k £ [| V (I - *0" * »L (*■.* + T^) i 

cos 8 <Pi + (■§ - <Pi) sin 8 cp x - ^ sin 2 qp x 

~ *L *S,k {" Bln 2c Pi " 4c Pi (i " «Pi)}] (57) 

71th radial loading, however, when the rings are not at- 
tached to the skin, we have: 

+ *L S | {«Pa cos 8 <Pi + (! - «Pi) sin 8 qp! - | sin 2qp x j- 

~ X l {tt sin 2^ - 4qp x (g - cp^j-J (5 8) 

Tor the practical computation of these coefficients, 
there may he used the chart in figure 11, giving the aux- 

iliary functions K (<P,X) and Z (<P,X) through which 
*S k and K 's k may 1)6 ex P ressei in simple form. The de- 
pendence on the angle cp x may, with good approximation, 
he expressed hy the following formula: 

*S,k foi) * \ <1 - c° 8 *q>i ) K S,k (45 °> (59) 

where K S t k» according to (57) and (58), respectively, on 
substituting the numerical values, Is: 

K S,k (45°) = Xg fk ^2.0833 + 3.6818 X L a (*S,i: + "f*^) 

- 5.5370 X L X Sfh: 10~ a (60) 

"a ■* —————— 

K (cp f A), for example, is the value of *s,k according to 
(57) for Xj, = Xg^jj. = X and. = 0. 
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for the caso of tangential loading., and 

■ r - - - g r( . ,. .... . , 

«S,k < 45 °> = *S f k (*sU + —a*) 12.0833 + 3.6818 Xj, 8 

- 6.5 370 X L ] 10~ 8 (61) 
for the case of radial loading. 

These formulas for the «s,k vo-luos of the circular 
ring may ho usod approximately for other ring shapes, 
lor the "exeontricities 11 X s k and X L mean valueB aro 

used; for (p 1 is substituted the angle which dividos the 

quadrant in the sane ratio as the stringer divides the 
oircumf erontial quadrants of the actual shell (<Pi =^ 

— i . A numerical example will make this apnroxl~ 

2 ti + "b 8 

mote mothod clear. Por an elliptical ring with ratio of 

axes 1 t £ and portions of circumf orence h x = h a = 

2 the evaluation of Kg ^ according to the Simpson rulo, 

if the contrlhution of the normal forces is neglected, 
gives : 



*S f k 



'0.000797 for X L = Xg tk = 0.95 
0.002173 for Xj, = X &fk = 1 



V. 

while the corresponding valuos for a circular "bulkhead 
with cp x = ^ are 0.000808 and 0.002282. Tho error for 

this narrow olllpse thorofore amounts only to ahout 5 per- 
cent. 

If tho panels have tho same dimensions axially, tho 
suhsoript k is droppod in the formulas and the displace* 
ment coefficients are then given in the form (31) with 

a» yr 4T>i f_*a_J? s* . 4^ / Ii >f s£ , X 

*l * Tl i = "~u~ \*rfs' ^ -.t.^r^l. ^ (62) 

(s* = mean wall thickness) 
and i% ' according to (57) or (58). 
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Tor the case of loading by a "convexing" force group „ 
Z 0 at the ond Bulkhead 0, the load coefficients afo 

forned according to (18) or (32), (33) fron the displace- 
ment coefficients of the redundancies. In the cat^> of mo- 
ments M k distributively applied in the elementary nah- 

ner, the load coefficients from. (6), (20), and (48) are: 

k -° ai L k S fk s 2f! / M ^s ltk+1 B 2fk+1 ^J 

(63) 

2. Six-Stringer Shells 

By assuming . additional stringers at the diametral 
points of the vertical axis, we obtain a 6-stringer shell 
(figs* 3 and 12). In addition to force group 1, there are 
two other types of redundancies, of which only the simple 
symmetrical force group (denoted in fig. 3 "by fc) 

is of importance. This force group occurs,' for examplo, 
in "bonding about tho 2 axis. 

Figure 12 shows tho unit state = 1. Thero is 

again a linearly decreasing force distrlhutlon in tho 
stringers. On account of tho symmetry about the y axis, 
there follows diroctly from the equilibrium at the stringer 
0: 

*l.k " - 2a^ ; = + 2^ (64) 

and further, from (4): 

The determination of the bulkhead stress requires in 
this case a simple, statically indeterminate computation. 
For a circa lar cylindrical shell with uniform stringer 

spacing (cp x = and with equal distances of all string- 

er centers of gravity from the center, this computation 
has been carried out elsewhere. (See reference 5, p. 464.) 
The bonding moments and the normal forces in the k* 11 hulk- 
head produced by the shears of the k* 11 bay, are: 
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for 0 < <p < ""'J 



(66) 



for * < <P < § 



H k = - X L (sin q> - | cob <p) for 0 < q> < ^ ^ 

Hfc " " 2a^ * L I C0B * f0r 5 - * - 1 



(67) 



J 



In a similar manner may be computed tho stress condition 
of rings not attached to the skin if radial load according 
to (27) is assumed. 

For tho displacement coefficients in the form (l7a,b, 
c), the w values for stringers and shoot according to 
(64), (65) are immediately obtained! 

Wt . + (ISit 1 ] (68) 

*B,k ■ = [a— + (r - O 8 b-^ 1 ] < 69 > 

' k L l f k 7 * S 2 f k G J 

Tho value of K s v ]c ln tlie °ase of genoral shape of cross 

section must again be determined by a numerical process* 
For a circular ring and with uniform stringer spacing 

^cp 1 = j\ tho value of «s,k may De computed on tho ba- 
sis of formulas (66) ,.(67) J. 

*S,k B *S,k [46.2136 + 57.4391 X L 8 (*s >k + % -fjS S ) 



- 103.0099 X L ^s.k] 10 ~* 



(70) 



As in the case of the "convexing" force group, this 
vnluo may also be used approximately for cross sections 
which deviate from the circular shape but whore tho 
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stringer spacing is about uniform. In loading the bulk- 
heads there are, in comparison with the convexing force 
group, twice as nany zero positions of the "bending moment. 
The contribution of the rings to the displacement coeffi- 
cient thus "becomes smaller as compared with the contribu- 
tions from the stringers and shear sheets.. The force dis- 
tribution then depends very little on the stress or the 
stiffness of the transverse stiffener walls. As the fol- 
lowing .numerical examples show, for shell shapes of such 
dimonaions as occur approximately in rionocoque fuselages, 
the statically indeterminate computation for the simply 
symmotrical redundancies may be carried out approximately 
under the assumption of perfectly rigid transverse stiffen- 
ere. Thus there is avoided the inconvenient computation 
of the bulkhead coefficients and a corresponding simplifi- 
cation is gained in the solution of the elasticity equa- 
tions. 

TTith uniform dimensions of the bays, the displacement 
coefficients are formed out of the K and u) values. Tor 
the two redundancies 1^=1 and T k = 1 these are, for 

convenience, collected in table of formulas 4. The dis- 
placement coefficients of the redundancies which occur in 
torsion, and bending in -the case of a 6-stringer shell with 
doubly symmetrical cross section, can then be obtained di- 
rectly from (31). By substitution of the values in tables 
of formulas 1 and 2, there is obtained the distribution 
due to a force group X 0 ani T Q , respectively, applied 

at the end bulkhead. 

The same redundant force groups hold also for an 8- 
stringer shell, obtained by adding two additional string- 
ers at the diametral points of the horizontal axis. In 
addition, there are then obtained two doubly symmetrical 
force groups which for the loading case considered are 
without significance, and one force group simply symmetri- 
cal about the z-axis which is effective in bending about 
the y-axis,. 

3. Shells with '.lore Stringers 

For shells with more stringers with double symmetry 
of cross section, the redundancies must first be so deter- 
mined that they have small effect on one another and their 
values may be determined with sufficient accuracy from the 
5-monber- principal equations. In simple cases this may bo 
done starting from the corresponding force groups with 
cyclical symmetry. 
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Ab an example, we consider a 12-st ringer shell with 
equal pane Is The ..force, gr.pups set up according to figure 
4, "break up into three "anti symmetrical , n two" doubly sym- 
metrical, and four simply symmetrical groups, of which 
only the antisynmetrical and those simply symmetrical 
about the y-axis need "be considered for the loading condi- 
tion investigated (fig. 13a,b). 

It is assumed that b x = h 4 ; b x = b a ; b B => b 3 ; 
Ba a s 3 ; and ff x = F 3 as, for example, is approximately 
the case with mono do quo fuselages. If the effect of tho 
bulkhead deformations - which with forco groups of highor 
ordor is small compared to tho offoct of the deformations 
of tho 3hoot and stringors - is nogloctod, then on account 
of the symmetry with respect to stringer 2 the antisynnet- 
rical force groups X^ k and X$ ^ are orthogonal to 
X 2> ]-» T he unknown forces c and d in X ljlc a 11 * ■ E 3 i k 
are now 30 determined that their mixed displacement coef- 
ficient contributions from stringers and sheet separately 
vanish. Then according to (31), any effect of different 
ttypos of force groups occurs /only thrqugh.the bulkhead de- 
formations, and therefore wi 1 * 3 ' = w^ x , 3 ' = 0. According 

to (12) and (14), these conditions are: 



1 



. 5jl ( 2o + i) (2d - i) 



(71) 



Prom tho so is obtained a quadratic equation for c- or d, 
the two roots giving, except for a factor, the same force 
groups Z]_ jj. and X.J ^ but in different order. For the 

special further case of b x =t T) a = b 3 = b- 4 = y|r, there is 
obtained for o and d: 

[ J*T~^ - ( ». - *i >] ] 

= g~ [ + i , i 8 + ct b - * )] 



21. 



(72) 



For the two force groups T 2 jj symmetrical 

about the vertical axis (fig. 13b3, there are determined 
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in similar manner the unknown forces c and d fron the 
conditions that the nixed displacement contributions from 
stringers and sheet vanish separately* The effect through 
the bulkhead deformations is then in both cases so small, 
as shown by numerical examples, that thoir values may be 
computed from the uncombined principal equations. 

The displacement coefficients from stringors and skin 
may for those force groups bo simply given by the general 
formulas (12), (14), (I7a,b,c). Whereas with groups of 
higher order the aspumption of rigid bulkhead may in gener- 
al be considered approximately true, in the case of the an- 
ti symmetrical force group ^ the bulkhead coefficient 

must bo taken into account. For ring bulkheads this coef- 
ficient can bo determined approximately from the Kg val- 
ues for the 4-stringer circular cylinder. The force group 
*l,k la composed of three simple convexing force groups 
at stringers 1, 2, 3. For these the Kg values are com- 
puted from, the formulas for the circa lar ring bulkhead. 
On account of b x = b 4 the values for the groups at string- 
ers 1 and 3 are equal to each othor, so that only the two 
magnitudes 

, \ t \ 4bi 4 (b x +b a ) 
» s (cp 1 ) and Kg (qp a ) for cp x = "y i qp a = ^ 

need be takon from (57), (58), or figure 11. TTe thon have 
approximately for the Kg v.".luo of the force group ^i^: 

*S = K S to a ) + 3c a Kg (tp x ) + 4c/its" (qp x ) *g(<P a ) (73) 

The accuracy will be checked with the aid of an exam- 
ple. For 7\g = 0.938, 7\ L = 0.9625, t3. = = g, 

<Pa = 3» there is obtained from figure 11: Kg (<p a ) = 
8.52 X 10*"*, K s (cp x ) =4.29 x 10""*. The approximate for- 
mula gives for c = 0.236: Kg k 14.94 x 10""*, while the 
numerical integration according 'to the general formula (16) 
givos tho moro accurate value 15.17 x 10"**. For practi- 
cal computation the error is insignificant. 

For shells with many stringers and with few axes of 
symmotry the setting up of forco groups which are approxi- 
mately orthogonal is vory tedious. A fundamental method 
is given in the Appendix. By combining several longitudi- 
nal stiffoners into one fictitious stringer, it is goner- 
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ally- possible to obtain the simple shell forms considered 
and thus a general .view of the ^03.0.0, .distribution "In the 
large." ' The approximate computation may likewise he re* 
stricted to small ranges of the system in tho neighborhood 
of "disturbance positions" (discontinuitios in loading or 
dimonsions, cutaways, etc.), in which oase arbitrary lin- 
early indopondont force groups are chosen as tho static, 
redundancies and the complete elasticity equations for 
these solved as in section III, 3. 



71. JTUHEHIOAL EXAMPLES AND COMPAHISOH TTITH TEST EESULTS 
l a Eorco Distribution .in a Six-Stringer Shell 
Effect of the Various Stiffnosses 



Tho effoct of tho various stiffnosses on tho forco 
distribution will bo invostigatod with tho aid of a sim- 
ple oxamplo of a 6-stringor circular cylindrical shell of 
equal bays with uniform stringer distribution, and tho 
roliabllity of approximate computations chockod. The bulk- 
head rings are assumed to bo attachod to tho skin. 



r = 40 cn; rg = 
E Q = 3.0 omS; E 1 a 
s x = s a = 0.06 cm; 

From table of formulas 4, tho id values, which are indo- 
pondont of the redundancies, aro: 

ujl = 0.333; (JUu = 1.088 x 10*; U)g = 58.8 X 10« 
Purthor: 

X S ■ ? = 0.925; At a ~ = 0.9625 

a) Applica t ion of a c o nvoxing forco group I Q = 1 
(fig. 14).— According to table of formulas 4 and figuro 111 

K L = 288; K B = 0.250; K S = 0.717 X 10~ 3 



Dimensions 

37 cm; r L = 38.5 cm; a = 36 cm 
4.5 cm 8 

J s = 3.0 on*- ^ « ^ 
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so that, k l u» L = 0.96 x 10 s ; K B cuj = 27.2 x 10 B 

Kg U) = 421.6 X 10 B 

Prop thoeo thore aro- obtained, according to tho formulas 
In tablo 2, tho values: A = 1.147; B = 0.341. ffo there- 
fore havo tho case (A - 3) < 1, with tho upper sign, 
ffo thon havo : 

\|f = | cosh -1 1.488 = 0.476, X = | (cos™ 1 ) 0.806 = 0.317 

The falling off to zero of tho convoxing force group X Q = 

1 for an infinitoly long sholl is thon ohtained for clas- 
tic and rigid end hulkhoads, roopoctivoly , from tho equa- 
tions : 

X k = e "*°'* 76k (cos 0.317 k + 1.055 sin 0.317 k) 

and 

Xfc = e ~°' 476lc (cos 0.317 k + 0.390 sin 0.317 k) 

For tho Uniting oace of perfectly rigid end hulk- 
heads : 

Iv = e--<P with cp = cosh" 1 27,g+B ' 84 = cosh -1 1.227 = 0.662 
* 27.2-1.92 

Figuro 14 shows the nunorical valuos of plotted 

against the shell length, tho points Doing joined hy 
straight linos. According to tho sholl nodol usod the 
forco distribution in one of the four loadod stringers is 
shown. For tho case of olastlc end "bulkheads the valuos of 
Xjc aro glvon for a 6-hay sholl. Largor deviations fron 
the valuos for tho Infinitely long sholl occur only for 
snail ond forces, so that oholls whoso longths aro approz- 
inatoly equal to tho circunf oroncoa , nay "bo conputed fron 
tho sinplor formulas for tho Infinitoly long sholls. 

The effect of the converting force group will ho snail— 
or tho stiffor tho transverso stiffenors conparod to tho 
stringors and shoot. If the stringers, howovor, aro vory 
stiff (k l (U l « o), then tho "bulkhead effect vanishos; 

l.o. , with finite longth sholls tho forces in this caso 
docroaso linoarly ovor tho ontiro longth. 

A further Uniting caso is that of a shoot rigid In 
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shear (k b " U) b « .0) and tho force distribution la shown in 

figure" 14 (thin lines') . ~ The" effect of the - shear deforma- 
tions on the force distribution appears to he groator the 
stiffer the "bulkheads; a great stiffness in shear gives a 
more rapid rate of docrease in the forces. The continuous 
distribution of the far ces computed according to table of 
formulas 3 gives at the transverse stiffenors about the 
sane values as obtained according to tablo of formulas 2, 
and therefore has not boon specially indicated. 

3?pr the statically indotorminato computation, it 1b 
important to know the effect of slight variations in the 
stiffnossos since the computation is made with estimated 
cross sections and skin thicknesses. Piguro 16 shows the 
forco distribution in a long sholl with rigid and ring 
whoro tho dimensions of bulkhoads, stringers, and skirt 
havo been variod individually . from tho main assumption a). 
It nay bo soon that tho force distribution is not very 
sonsitive to those changos. From cases c) and d) it nay 
bo concluded that tho total distribution after buckling of 
tho skin does not chango nuch, since tho of feet s duo to 
the anallcr contribution of tho skin in supporting tho 
longitudinal forces and tho doeroaso in tho shonr stiff- 
ness approximately offoot oa.ch othor. Upon this force 
distribution there 1b still to be superposed, however, the 
states of stress due to the tension fields themselves. 

' b) Applicatio n of a synnetrlcal forc e g roup T 0 = 1 
(fig. 16). The type of loading shown in figure. la by a 
bending force group may, according to section II, 3 be 
reduced to a simply synnetrloal force group T 0 . Pron ta- 
ble of formulas 4, we find: 

KL = 360, K B = 0.1465, Kg = 1.085 X 10~" 8 

fron which A = 7.50, B = 5.77. We thus have the case 
A - B > 1 with upper sign, and the falling off to aero, 
for an infinitely long shell is aperiodic, for example, in 
the case of elastic end bulkheads ,. according to the formula: 

Tjj. = a"* 1 ' 6381 * (cosh 0.573 k + 1.812 sinh 0.573 k) 

The results of computations for the different limiting 
cases of bulkhead stiffness are plotted in figure 16 against 
the shell length and the points joined by straight lines, 
the force distribution in stringer 0 being represented. 
There is also shown the limiting case (Xl = >>s ~ 0 ) >for 



50 



1T.A.C.A. Technical Memorandum Ho. 8661 



which Kq = 7.35 X 10" 5 .. This case may also he interpret- 
ed with respect to the original dimensions as a reduction 
in the bulkhead stiffness by. about 1/7. There has further 
been computed the force distribution in the limiting case 
of sheet rigid against shear. 

Since the force group T 0 is of "higher order" (see 
II, 2) than the convoking group X 0 , it reduced to zoro 
over a shorter sholl length. The force distribution do- 
pends only on the stiffness of the transverse stiffoners, 
particularly if there is a strong end bulkhead at tho load- 
ing end. The dependence on the shear stiffness is consid- 
erably greater, however, than in tho case of the convexing 
force group. As an approximating assumption for simplify- 
ing the computation, it is accordingly permissible to neg- 
lect tho bulkhead deformations (*s = 0) * n t ^ 10 static- 
ally indoterminate computation for the simply symmc trical 
force groups. 



2. Computation of a Test Sholl 

The computation procedure developed trill be applied 
to a circular cylindrical 18-stringer shell, whose states 
of stress under an applied bending and convexing force 
group at 4 stringers, has been determined oxporimentally , 
as presented in detail in a simultaneously appearing papor 
by S. Schapitz and G. Krumling. (Soo reforenco 1.) The 
description of the shells and tho method of conducting tho 
toots may thoreforo here bo disponsod with. We furthor 
limit oursolvos to tho computation of the longitudinal 
stresses in tho stringors bofore budcling of tho skin, 
takes place, and comparo them with tho measured values. 

a ) Applicat io n of a bending force group .- Tho actual 
arrongoaont of tho stringors is shown in figuro 17a for 
one quadrant. At strin^or s [7] , [12], [3], [16] 10 , a puro 
bending momont about tho transvorso axis is applied by 
four concentrated forces P. T7o consider first tho closod 
shell without the cutaway portion at tho loading sido and 
idoalize the systom to a 12-stringor sholl by combining 
stringors [10] and [ll] into one stringor 1 of doublo 
the cross— sectional area. Tor tho skin thickness thoro is 
assumod a uniform valuo of s = 0,06 cm and half tho skin 



Tho stringor and bulkhead notations taken from the papor 
by E. "Schapitz and G. Krumling aro donotod by brackots. 
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atrip of the neighboring sheet panels is added as offec- 
tivo width to each stringer soction. Thqsq_..ef f Qctiyo 
•stringer sections are shown in "figure . 17a. We' then have: 

4 (P x 7 X B + I a 7a 8 + * s 7 3 8 ) " 2 - 50 x 1q4 0m4 

equal with sufficient approximation to the moment of iner- 
tia with respect to the . transverse axis of tho full shell 
in tho last bayB. Tho applied bending force group is split 
up as ahbwn in figure 5. The component forces of the force 
group corresponding to the linear stress distribution: over 
the cross section and having oqual moment are: 

P,( q) = 0.383 P; Pj°> = 0.389 P; P^ 0 ^ = 0.097 P 

As redundancies there occur only the simply symmetrical 
force groups T i t k and T 2 indicated in figure 13b. 
From the conditions for the vanishing of tho mixed displace- 
ment coefficients from stringors and shoot there are ob- 
tained the values for tho componont forces c, d, c, d 
shown in figure 17b (solution of a quadratic equation). 
The characteristic forco group shown in figuro 5 as tho 
difforonco botweon tho actual and linearly distributed mo- 
ment is combined from tho redundant longitudinal forco 
groups in tho following nannor: 

Characteristic forco group = 0.431 P (T 1>0 = l) 

- 0.262 P (T a>0 = 1) 

Tho manner of decrease of T x 0 = 1 and T a 0 = 1 is com- 
puted under tho approximating assumption of rigid bulk- 
heads - which assumption is permissible according to tho 
numoricr.1 oxample previously given. The bulkhead spacing 
is a = 36 cm and with tho valuos determined according 
to (28) and. (30): 

K L = 2 ? 5 3 x 10 B ; K B U) B = 32,9 X 10 a for T ljk = 1 
K L idl = 2.86 X 10 a ; K B (Uj = 12.84 x 10 a for T 2>fc = 1 

thorp is. obtained from table of formulas 2: 

T l f k - e- 1 ' 001 . k; T 2fk = 0 - 1 ' 898 k 
Prom these values are obtained the stringer forces due to 
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the characteristic force group^ and "by superposition with 

the above values of P^ 0 \ P B ° , P 3 ° , there is obtained 
the final force distribution plotted in figure 18. 

In order to take into account the effect of the cut- 
away portion, this incomplete hay is separated from the _ 
full shell and a further simply symmetrical force group T 
is assumed to act as static redundancy at stringers 2 and 
3. For the condition T = 0 there is obtained in the 
full shell the force distribution previously computed, in ■ 

ihe main spars of the cutaway portion the constant, force 
P. The stress condition 7 = 1 in the full shell is re- 
stricted ap£roximately to the first hay. For the determi- 
nation of T, the various skin thicknesses must he taken 
into account. Furthermore, the skin in the first hays 
does not contribute to the support of the longitudinal 
forces in the' same degree as was previously assumed in the 
case of the complete shell. After computing the displace- 
ment coefficients according to the general formulas in 
section III, 2, using the values indicated in figure 18 
for the dimensions (bulkheads assumed rigid), there is ob- 
tained ? = 0.0606 P, and from this the force distribu- 
tion given in figure 18. 

For comparison with the measured stresses the mean 
stringer stresses must be determined from this force dis- 
tribution. Since from bulkhead [d] on the stress dis- 
tribution obtained is approximately linear (see fig. 18), 
it is possible from this position on to use the values 
for tho stringer cross sections given in figure 17. At 
tho cutaway portion and at the first bays, the skin con- 
tributed only imperfectly to the support. (See the paper 
by Schapitz and Krumling, reference 1.) At the cutaway- 
portion tho mean cross-sectional areas given in figure 18 
apply; at bulkhead [f] in the full shell, wo have ap- 
proximately: 

F x = 1.95 cm 8 ; F s = 3.10 cm a ; F 3 = 1.30 cm 8 

From the Be values up to those of figure 17a at bulk- 
head [d] a linear rate of increase is assumed. In fig- 
ure 19 are plotted the stresses computed with these cross 
sections from tho force distribution in figure 18, as well 
as the |( values given in figure 7 of tho work by Schapitz 
and Krumling, the values being divided by an initial 
stress of 474 kg per cma a t the loading side. The stress 
distribution is in satisfactory agreement with the stat- 
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ically lndptorminate Computation. The assumption of rigid 
hulkhoada in the computation loads to a more rapid state 
of docrease of the longitudinal strosses in the main spars 
and to a more rapid' rate of rise in the intermediate 
stringers. Through the combination of stringers [10] 
and [11] the deformations of these stringers occurring 
in tho first hays of the oomplete shell are not taken into 
account, thus resulting in a. greater stress in stringer 
[11] than in stringer [10]. Furthor details of the 
stress condition (shears, bulkhoad loading, etc.) are given 
in tho experimental report referred to. 

b) Application of a, oonvexlng foroe group.- For com- 
puting, the force distribution resulting from the applica- 
tion of a "convexing" forco group at " the four main spars, 
the test specimen is idealized into a 12-stringer shell In 
a different manner. Since the intermediate stringers in 
tho neighborhood of the upper and lower diametral points 
are hardly loaded this computation is restricted to the 
main spars and the adjacent intermediate stiffeners. and 
is "based on a somewhat differont circumf orential spacing, 
as shown in figure 20a. The offect of the cutaway portion 
is neglected. For an approximate computation, we consider 
first a system of equal bays with skin thickness s «= 0.06 
cm and a bulkhead moment of inertia Jg = 3 cm 4 (this is 

approximately the computed value for the section of the 
bulkhead rings [e] and [f] attached to the skin at the 
loading side. ) To the main spars and intermediate string- 
er sections (nee fig. 17a) there is added an effective 
skin strip of 12 and 14 cm, respectively, and there are 
thus obtained the effective stringer sections shown in 
figure 20. 

Of tho throe ant 1 symmetrical characteristic force 
groups at a 12-stringer shell (fig. 13a), only and 
"^Sjk' ^ UQ *° *ke convexing forco group, will be effec- 
tive. For tho component forces c and d, there are de- 
termined tho valuos indicated in figure 20b from the for- 
mula' (72) in section 7, 3. The mixed displacement coeffi- 
cients duo to stringers and sheet then vanish. For the 
applied convexing force group Zp, the following relation 
is true: 

Xp = 0.780 P (X 10 m 1) + 0.220 P U Zf0 = l) ' 

The forco distribution duo to Z x 0 = 1 and X 3j0 = 1 is 
detormined from independont dlfforonce equations, as shown 
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in the first numerical example. The bulkhead coefficient 
Kg for the redundancy ^ j. Is already contained in the 

numerical example for the approximate formula (73). Fig- 
ure 21 shows the results for different assumptions. The 
characteristic force group X 3 0 = 1 decreases more rapid- 
ly than Z, 0 = 1, since the former is of higher order* 
i 

It nay ho determined with sufficient accuracy according to 
tahle of formulas 2 under the assumption of Infinitely 
long shell and rigid transverse stiffoners. For the char- 
acteristic force group X-^ k these assumptions are no 

longer admissible; tahle of formulas 1 must thereforo "bo 
used and the effect of tho "bulkhead dof ormations taken in- 
to account. Thoro was, furthermore, investigated the of- 
fect of the two redundancies duo to the "bulkhead deforma- 
tions, "by computing the mixed displacement coefficients 
and solving the complete elasticity equations. This effect 
is very small, the redundancies "being practically orthog- 
onal • 

Slnco the variation of the redundancies Xi,k depends 
essentially on the "bulkhead stiffness, a more accurate com- 
putation was made for the actual dimensions "stepped" in 
the longitudinal ■ direction. Tho effective moment of in- 
ertia of tho "bulkheadB at the sholl attachment was ob- 
tained by deflection measurements for tho case of diamet- 
rically situated concentrated forces: 

Jg = 6.5 cm 4 for bulkhead [o] and [f ] 

J s = 0.9 cm* for bulkhead [a] to [d] 

For the bulkheads [o] and [f] attached to the skin, tho 
bulkhead coefficient undor tho assumption of tangential 
load is: 

Kg = 14.94 X 10"" 4 

while for bulkheads [a] to [d] under the assumption of 
radial load (according to fig. 11 and approximate formu- 
la (73)): 

(t g ' = 23.9 X 10"* 4 

Furthermore, there was taken into account the variable 
skin thicknesses (0.08 cm in bay 1, 0.06 cm in bays 2 to 
4). Under these assumptions there are obtained, according 
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to the general procedure in section I'll, the values for 
Z l jfc' Plotted as - 8 heavy- line in figure. 21. In spite of 
the weakor .bulkheads at the ond of the sholl thoro is a 
more rapid rato of decrea'so than undor tho assumption of 
equal hoys. This is due to the essentially smaller stiff- 
ness of the system from hays 2 to 4 as oompared with the 
first "bay. 

With those more accurate values for ^i^k a 314 * n0 
approximate values of X3 ^ for rigid bulkheadSi the forco 
distribution shown in figure 22 was computed. lor compar- 
ison with the experimental results thoro wore further de- 
torminod the stresses duo to this force distribution. Tor 
this purpose there was assumed a linear rate of inoroaso 
botwoon tho estimated reduced areas ?i = 1«35 cm 8 , 3? B = 
3.10 on? at bulkhead [f] &ad tho stringer cross so ct ions 
from bulkhead [ d] on, according to figuro 20. Since tho 
offoct of the cutaway portion is neglectod, the comparison 
can only be made with the measured stress values in the 
main spars and in the intermediate stringers [11] and [173 
(fig. 22). Tho experimentally detorninod strossos in. tho 
main spars decroaso somewhat more rapidly than according 
to the computation, ospecially from tho second bay on* Ao 
a rosult, the computed shears (nee fig. 20 of tho papor by 
Schapitz and Krunling) in tho first bay aro snaller, and 
in tho laot bays aro groator than tho oxporinentally do- 
terninod values. Tho reason for this lies in the nogloct- 
ed bonding stiffness of the nain spars, as a result of 
which tho bulkhoads at tho ond of tho sholl rocoive addi- 
tional strossos in tho samo dirootlon as those from the 
shears. Thoso bulkhoads therefore act still moro weakly 
than was assumed in tho computation, and the roduction is 
concentrated noro strongly on the stiffer first bay. In 
spito of tho neglect of these factors, the method devel- 
oped yiolds sufficiently acourato values. 

Further conclusions as to the state of -stress need 
not here be considered as those will be found in tho work 
of Schapitz and Krumling, already roferrod to. Tho effoct 
of modifiod teat conditions may bo roughly estimated from 
the investigations in the first numerical example. In 
particular', tho tests confirm -the conclusion' drawn from 
figure 15, that the force distribution as a whole under- 
goes only a Blight chango aft or buckling of tho skin. 



56 



1J.A.C.A. Technical Momorandun Ho-. 866 



7I.I . SUMMARY . . 



Longitudinally and transversely stiffened shells are, 
under the assumption of constant shoar flow' (shear times 
thickness) in each of the sheet panels, treated as systems 
with finite static redundancies* The procedure of the 
statically Indeterminate computation is given for cylin- 
drically shaped shells under conditions of loading "by con- 
centrated forces, moments, and transverse forces. As stat- 
ic redundancies longitudinal force groups are introduced 
at the intermediate transverse stiffeners and restraints 
betwoon tho hays, theso force groups "being so chosen' that 
only groups of the same kind affect each other longitudir 
nally — as a result of which the elasticity equations "break 
up into independent 5~mom"ber partial systems. The manner 
of sotting up of these "orthogonal" characteristic force 
groups is given for several simple shell shapes* Such 
force groups are In general possihle only for a sufficient- 
ly large number of axes of symmetry, "but otherwise only, if 
tho mutual offect through, the doformations of tho trans- 
verse walls or tho shear sheets or also of the stringors 
is negloctod. (Soe VIII, Appondix. ) It is shown, with 
the aid of oxamplos, that neglecting of the transvorso 
stiffonor doformations Is admissible for tho setting up of 
orthogonal force groups. The further result is obtained 
that the transverso stiffoner deformations are of impor- 
tance only for tho distribution of the force groups which 
r.ro in equilibrium over tho entire shell circumference, 
(for oxample, anti symmotrical force groups in torsion)* 
For force groups of "higher" ordor - such, for examplo, as 
occur in bonding - tho transverse walls may be assumed as 
rigid. 

In tho case of equal bays, the elasticity equations 
arc solvablo in finite form. For tho loading of the shell 
by concontratod forces, these solutions aro given In ta- 
bles of formulas. Similarly, thore havo beon colloctod 
in one tablo, tho displacement coefficients for 6-stringer 
shells of f qual bays. TTith tho aid of those tablos tho 
most important disturbances from the olemontary force dis- 
tribution, such as occur on application of concentrated 
forces at cutaway portions, may bo approximately deter- 
mined for the usual sholl shapes. 

The practical computation procedure is clarified with 
the aid of a simple numerical example, and the effect of 
tho stiffnesses on the force distribution, investigated. 
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Further, there is computed according to the proceduro pro- 
poaod, the stringer stresses in a test shell with axial 
f arceo' applied at four points and tho" stresses - compared 
with the experimentally determined values. Satisfactory 
agreemont is obtainod "between the computed and measured 
values. Tho same procedure, using the corresponding dis- 
placement coefficients, may he applied to othor multi- 
stringer systems (multispar frame wings, airship hulls). 



VIII. APPENDIX 
Orthogonal Poroe Groups for Cylindrical Shells 
with Arbitrary Cross Sections 



Tho sotting up of uniform orthogonal characteristic 
force groups for shells of equal hays with arbitrary non- 
symmetrical cross section and many stringers, will ho moro 
closely investigated. If, in tho mixed displacement coef- 
ficients of tho forco groups at tho same transverse stif- 
fenors., tho contributions from the stringers, sheet, and 
bulkhoads vanish separately, thon according to (31) all. 
mixod displacement coofficionts of force groups of differ- 
ent kinds aro equal to zero and conversely. It is thoro- 
foro sufficient, for tho sotting up of orthogonal groups, 
to considor a system- of two bays and at the intermediate 
transvorso stiffener, to doternine tho redundancies X x , 
X St ■ ■■•! In-3 is such a aannor that all mixod displace- 
ment contributions 6,^ L ,| , 65. B ,? , and 8,^ S ,) vanish. Thore 

aro 3f n 2 3 J = | (m - 3) (m - 4) different portions of 
this kind. 

Starting from arbitrary linoarly independent force 
groups' Zj, , X a X n _ 3 and forming from those by a 

linear transf orma'tion with determinant different from zoro , 
tho requirod forco groups X a X a-3 : 

_ m-3 

Z f* = W I X ^l?. Cpi = 1, 2, .... n - 3) (74) 

thore aro only (m-3) a - (m - 3) = (m-3) (m-4) os- 
sontial constants c^ v available since oaoh multiplo of 
a forco group may be oho son as a unit state. The -| (m - 3) 
(n — 4) oquations of condition for tho vanishing of 
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3 (L) 8^ 3 ,^,'and 6^ S ^ cannot all, therefore, as a rulo, "bo 

satisfied. In special cases - as, for example, with cyc- 
lic synnotry - those equations nay ho independent of ono 
another, so that the froo valuoo of tho ahove transforma- 
tion aro sufficient. In gonoral, howovor, the equations 
are independent, and the orthogonality condition of the 
force groups cannot he obtained for arbitrary nonsymmetri- 
cal cross soctions. 

If the offect of ono of tho throe stiffnesses of 
stringors, shoot, or bulkheads on tho coupling of tho dif- 
ferent redundancies is noglectod, then tho (n - 3) (n - 4) 
osscntial valuos &\i t v of tho transformation aro suffi- 
cient for satisfying tho ronaining 2 ( n g 3 ) = (m - 3) (m - 4) 

conditions. In tho work of Wagner and Sinon (reference 3), 
the offect of the 3hoar stiffness is not takon into account, 
so that tho determination of orthogonal independent stress 
groups is possihlo as characteristic solutions of a linoar 
integral equation with symnotrical nuclouc. As is ohown 
by tho numerical exanplos for stiffoned sholls, the offoet 
of tho bulkhoad def ornations on tho force groups of higher 
order is small, so that it is convenient to noglect tho 
offoct of tho bulkheads - i.o., to dotermino the force 
groups Xjj, in such a manner that tho mixod coefficients 
(L) (b) 

8., „ and fi M vanish. The strain energy of tho systom 

conposod of two panols nay bo exprossod in terns of tho 
redundancies X\ , .... Xri -3 with tho corresponding dis- 
placonont coefficients 6n,v ^ke quadratic portion in X 
consists, on noglecting the bulkhoad deformations, of tho 
two positive definite quadratic forns: 

H (L) m 1 E ? (L) s , (B) m 1 s ? (B) j ^ 

whoro the X's now donote nunbors, nanoly, nultiplos of 
the unit statos of the redundancies X. According to a 
theorem on quadratic forns (roferonco 14), thoro thon al- 
ways exists a linear transformation: 

by which and H^ 3 ^ aro transformed into oxpros- 
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slona with purel7 quadratic terms: 



' s CiO = i V ,(Lj ,7' h (b) . i -j» A3) x . (re) 
o 

She coefficients ^ of the transformation are the 

(m - 3) particular BOlutione of the homogeneous systom of 
equations: 

HI (Z^l - Pli 8$) 0,,^ - 0 (X - I, E m - 3) (76) 

T7hoso determinant, set equal to zero (an equation of the 

(m - 3) 4 dogroo), has the characteristic numhors 

(p. = 1 , 2, . . . f n - 3) ao roots. If tho transformation 

m- 3 _ 

is nor applied to tho unit stato of tho force groups, tho 
nixod displacement coofficionts from stringers and shoot 
for the fcrco groups I vanish, since for 6^^, ne have 
the transformation: 

8 ».» ■ pil. vi <77) 

and for H^ L ^ wo have, aftor substituting from (74') and 
(77): 

« (t> - * A ^ (S (f x.)- 

" I ^ *» [p?4 •»>»] 
"•5 ^ 8 $ ^ * - I 5 < L ' V 

so that, 

(L) (L) (L) - _ _# 

8 H,u " d n fo * 14 " v ' «|*.-» = 0 for ^ £l> 



Similarly, for 8 



(B) 
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Tho dotornination of orthogonal characteristic func- 
tions by noglocting one of tho throo essontial stiffnesses 
with rospoct to thoir mutual coupling effects is thuo ro- 
duced to the problem of the simultaneous tranaf ornation of 
two positive dofinite quadratic f orris into expressions 
with only square t orris. In tho case of n stringers, 
this requires tho complete solution of a characteristic 
equation of the (n - S)& degree and tho computation of 
the corresponding particular solutions of the linear homo- 
gonoous system of equations (76). These algebraic oper- 
ations .for a finite numbor of variables correspond to tho 
solution of the linear homogeneous Integral equation in 
tho work of Wagner and Simon. 



Translation by S. Eeiss, 
National Advisory Committee 
for Aeronautics. 
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figs. 1,2,3 




Figure 2.- 

Sketch of 
•hell 

indicating 
notation 
-used. 




k+f n 
Bulkhead number 

a. Concentrating 
bending forces. 




b, Concentrated 
distorting 
forces. 





Figure 1.- Examples of disturbances 

Force group k from the elementary 

» stress states. 



Force group Xg^ 
Force group 

Figure 3.- Force groups for 
six-stringer 
shell. 
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Tigs. 4,5,6,8 



r 1 



= 2 






r/w-Qcpsp—fr^^^ group + characteristic force 

group. 






Bending force group = linearly distri- 









rift 
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P<°> = P 



F, 1A j/» 
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Simply-symmetrical redundancies 
Figure 5.- Example of the splitting 

up of applied longitudinal 
concentrated forces. 




Figure 4.- Force groups for 12-etringer 
shell with cyclic symmetry. 





Figure 6.- Stress 

distribution 
in principal bay 
system due to X^jj. ■ 1. 



Figure 8.~ Loading of "bulkhead by radial forces 
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Tigure 7.- Scheme for elasticity equations for n = fi and m = 5. 
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Tigure 11.- Determination of bulkhead coefficient's^ 

applied convexing force .group Xjp 5= I, 



for the 



for tangential loading (bulkhead) 
for radial bulkhead loading 
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Figs. 9, 10, 12, 13 




t, 

Bu/khead A 



Figure 9.- Stress condition due to applied convexing 
force group X k = 1. 




a, Symmetrical with respect 
to center. 
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-pi 















ft ' s, 
b. Simply symmetrical 

Figure 13.- Characteristic force 

groups for 12-stringer 
shell with doubly-symmetrical 
cross section. 





Figure 10.- Section of a 

circular 
cylindrical shell. 




Neutral' t] 
ring axis Bulkhead k 
Figure 12.- Stress condition due to simply-symmetrical 
characteristic force group Y k = 1. 
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Figs. 14,15,16,1? 



Infinite- 
ly long < 
shell 



All bulkheads 
elastic. 



b, End bulkhead 

rigid. 

_ c, All bulkheads 

rigid. 
Thin liner: limiting 
case of sheet rigid in 
shear. 




Bulkhead 
number k 

Figure 14.- Falling off of applied 

convexing force group 
Xo = 1 over shell length. 
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Figure 15.- Effect of the various 

stiffnesses on the 
distribution of the applied force 
group X 0 = 1. 




Figure 16.- Falling off of simply- 
symmetrical character- 
istic force group T 0 = 1 over 
shell length. 




b, Effective redundant force groups 

Figure 17.- Bending force groups 
applied to the 

test shell. 



N.A.C.A. Technical Memorandum No, 866 



Figs. 18,19,20,21,22 



Ttesf values 



Simplified 




it 



._L_ 



Figure 18.- Force distribution in 

the simplified shell 
due to the bending force group 
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Figure 19.- Comparison of computed 

and measured stringer 
stresses due to the bending force 
group. 



Solution for equal bays .from 
elasticity; equations 

l -> T Bulkheads 

Independent ^ -/elastic 



difference 
equations 



a, Simplified shell shape 






Rigid 



*i t k for stepped 
dimensions 



b, Effective force groups 
Figure 20.- Convening force group 
applied to the test 
shell. 



Figure 22.- Force distribution 

and comparison of 
computed and measured stringer 
stresses due to convexing 
force group. 
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Figure 21.- Distribution of the 

antisymmetrical 
characteristic force groups 

^ - 1 and X 3 k = 1 in the 
simplified shell. 
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